UNIVERSIDADE FEDERAL DO RIO GRANDE DO SUL 

INSTITUTO DE FÍSICA 

DEPARTAMENTO DE ASTRONOMIA 


Estudo teórico dos 
perfis de linhas do Hidrogênio 

perturbado por colisão com prótons 


Márcio Gabriel dos Santos 


Tese de Doutorado realizada sob ori¬ 
entação do Prof. Dr. Kepler de Souza 
Oliveira Filho e apresentada ao curso de 
Pós-Graduação em Física da Universidade 
Federal do Rio Grande do Sul, como req¬ 
uisito parcial para a obtenção do grau de 
Doutor em Ciências. 


Porto Alegre, RS, Brasil 
Outubro de 2012. 



Ex umbris ad lucem! 



Agradecimentos 


Quero agradecer primeiramente ao meu orientador professor Dr. Kepler de Souza 
Oliveira Filho, por esta longa orientação, pela sua paciência, seu apoio, mas princi¬ 
palmente pelo incentivo e confiança. 

Agradeço em especial àquela que sempre me apoiou incondicionalmente, que 
apostou em mim mais do que ninguém e que seguramente é a que mais compartilha 
da minha alegria: Elizandra Martinazzi. 

Por fim, agradeço aqueles que, de alguma maneira, contribuíram com este tra¬ 
balho e que me incentivaram muito! 



III 


Resumo 

Estudo teórico dos perfis de linhas de Hidrogênio 
perturbado por colisão com prótons 

O foco deste trabalho é o estudo da interação do átomo de hidrogênio com 
prótons, criando momentaneamente a molécula H^, determinando os potenciais do 

para níveis atômicos mais altos (n>3), bem como analisando todas as possíveis 
transições que podem ocorrer. O cálculo do momento de dipolo elétrico do 
para todas as possíveis transições que encontramos é de suma importância para 
determinar as linhas espectrais observáveis. A aplicação de tais dados teóricos terá 
sua inserção no estudo dos perfis de linhas de estrelas anãs brancas DA, com o 
intuito de obter melhores ajustes dos modelos de atmosferas estelares. Para isso 
calculamos todos os potenciais relevantes para as transições das séries de Lyman e 
Balmer até n=10 (número quântico principal) do Hj e os respectivos momentos de 
dipolo elétricos para tais séries. Na literatura encontra-se somente os valores para 
Ly a, Ly/3, Lyy e Ha, mas mesmo estes, com aproximações mais grosseiras do que 
as que calculamos. 

As leis da física permitem o cálculo de um modelo de atmosfera estelar num dado 
instante de sua evolução e o confronto com os dados observacionais. Entretanto, por 
melhor que seja o modelo, ainda não se obteve um resultado plenamente satisfatório 
entre teoria e observação. Ainda existem muitas incertezas nos modelos estelares, 
principalmente nos cálculos de opacidade e convecção, quando comparamos aos da¬ 
dos observacionais, o que implica na necessidade de cálculo mais detalhado da física 
que rege tais modelos. 

O objetivo principal desse trabalho é calcular os potenciais de interação das 
partículas constituintes da atmosfera e analisar as linhas espectrais das estrelas que 
já estão na última fase de evolução, as anãs brancas com uma atmosfera pura de 
hidrogênio, chamadas de DAs. 

Um dos melhores métodos para a determinação da temperatura efetiva e obtenção 
da gravidade superficial de anãs brancas, é a comparação do perfil da linha obser¬ 
vado, a forma da linha espectral com o predito teoricamente. A eficiência desse 
método está no fato que o perfil da linha teórico é extremamente sensível à variação 
dos parâmetros atmosféricos - os perfis de Ha e H/3 são bons indicadores de temper¬ 
atura efetiva (T e y), enquanto os perfis das linhas mais altas são bons indicadores da 
gravidade superficial. 
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Entretanto, a técnica supostamente simples possui um grau elevado de difi¬ 
culdade, pois precisamos saber quais são todas as causas físicas responsáveis pela 
variação dos parâmetros atmosféricos que determinam o perfil da linha. Este é o 
objetivo de nosso trabalho, o cálculo do perfil de linha de H perturbado por colisões 
com prótons. 




V 


Abstract 

Theoretical study of the line profiles of the 
hydrogen perturbed by collisions with protons 

The subject of this work is the study of the interaction of the hydrogen atom with 
protons, temporarily creating from this interaction the H^ molecule. We estimated 
the potential for all leveis with n < 10 and analyze all possible transitions that 
can occur. We also calculated the electric dipole moment of H^~ for all possible 
transitions that we find, necessary for the line profile calculation. The application 
of such theoretical data will be in the context of the line profiles of DA white 
dwarfs, as refinement in the models of stcllar atmospheres. For this, I calculated all 
the potential relevant to the transitions of the Lyman and Balmer series to n=10 
(principal atomic quantum number), and their electric dipole moments. In the 
literature we find only Lya, Ly/3, Lyq and Ha calculations, but even these at lower 
resolutions and with the asymptotic approximations. 

The laws of physics allow the calculation of a stcllar atmosphere model for a 
particular moment in its evolution and the comparison with observational data. 
However, no matter how good the model is, it has not yet achieved a fully satis- 
factory result. There are still many uncertainties in the stcllar models in relation 
the observational data, showing the need for a more detailed analysis of the physics 
that governs such models. 

The main objective of this thesis is to analyze the spectral lines, calculating the 
potential of interacting particles constituents of the atmosphere of white dwarf stars 
with pure hydrogen atmosphere, called DAs. 

One of the best methods for analysis of spectral lines, determining the effective 
temperature and obtaining surface gravity of white dwarfs with hydrogen lines, is 
the comparison the observation of line profile, the sliape of the spectral line, with the 
theoretically predicted. The efficiency of this method is based on the fact that the 
theoretical line profile is extremely sensitive to variations of atmospheric parameters. 

However, the supposedly simple technique has a high degree of difficulty, because 
we need to know all the causes responsible for the determination of the physical 
parameters of the atmosphere affecting the line profile. This is the purpose of our 
work, the calculation of the line profile of H perturbed by collisions with protons. 
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Capítulo 1 
Introdução 


Ao imaginar que as linhas espectrais são formadas pelas transições eletrônicas entre 
níveis com energias muito bem definidas de átomos e moléculas, concluir-se-ia que 
essas linhas deveriam ser funções delta da frequência e apareceriam como linhas 
infinitesimalmente estreitas nos espectros estelares. Entretanto, muitos processos 
tendem a alargar essas linhas, de modo que elas desenvolvem uma forma -perfil- 
característica (e.g., Margenau & Watson, 1936). Os efeitos que alargam as linhas 
espectrais, intrínsecos aos átomos e extrínsecos, são o alargamento natural e as per¬ 
turbações produzidas pelas partículas vizinhas no gás, bem como pelo movimento 
dos átomos. Mas cada um desses efeitos produz seu perfil característico. As formas 
características são muitas; por exemplo, alguns perfis tem a linha central alargada e 
pequenas asas, outros possuem uma linha central estreita e grandes asas alargadas. 
A representação adequada desses efeitos permite o cálculo do perfil de linha obser¬ 
vado e, no processo, revela dados sobre as condições físicas nas camadas da estrela 
onde se origina o espectro, ou seja, comparando-se perfis calculados e observados, 
pode-se obter informações sobre abundâncias relativas de elementos químicos e seus 
estados de ionização e excitação, temperaturas efetivas e densidades das regiões de 
emissão do espectro. Portanto, para entender como as linhas se formam, é necessário 
conhecer os mecanismos de alargamento. A largura e a forma das linhas espectrais 
são de interesse em pesquisas físicas e astrofísicas porque: 

• fornecem informações sobre a temperatura, densidade e composição química 
existente na fonte do espectro; 

• o perfil de linha é importante na determinação das características físicas da 
atmosfera estelar, ou do plasma emitente para experimentos em laboratório, 

• a interação da radiação com átomos fornece informações para o cálculo do 


2 



Introdução 


3 


perfil de linha. 


1.1 Mecanismos de alargamento das linhas espec¬ 
trais 

A representação correta de um espectro implicará no conhecimento de todos os 
processos físicos que produzem esse espectro. Portando, para entender como as 
linhas se formam, é necessário conhecer os mecanismos de alargamento. Dentre os 
processos de alargamento das linhas espectrais, três são os mais significativos: 

1. alargamento natural; 

2. alargamento Doppler, 

3. alargamento colisional ou devido à pressão. 

O alargamento natural é devido ao tempo de vida finito dos estados envolvidos em 
uma transição eletrônica. É geralmente importante somente em baixas temperaturas 
e pressões, mas é o único mecanismo que persiste mesmo para um átomo estacionário 
isolado. 

O alargamento Doppler é devido a alteração das frequências da radiação ab¬ 
sorvida ou emitida por um átomo em movimento e, portanto, ocorre em um gás a 
uma temperatura não nula, pelo movimento térmico das partículas. 

E o alargamento colisional é devido às colisões de átomos com perturbadores que 
podem ser outros átomos, íons ou moléculas. 

1.1.1 Alargamento Natural 

A aproximação básica do alargamento natural está centrada no princípio da incerteza 
de Heisenberg. Ao considerar um átomo em um certo estado, o período de tempo 
At que ele pode permanecer nesse estado está relacionado à incerteza da energia 
A E desse estado por: 


AEAt « —. (1.1) 

27r 

onde h é a constante de Plank. 

Se houver um grande número de estados para os quais o átomo pode fazer uma 
transição, então a probabilidade de que isso ocorra é grande, logo o At será pequeno, 
e a incerteza da energia do nível será grande. Uma grande incerteza na energia de 
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um estado implica que um grande intervalo de frequências pode estar envolvido na 
transição associada a esse estado. Consequentemente, qualquer linha resultante da 
tal transição será alargada, conforme figura (1.1). 



Figura 1.1: Alargamento natural de acordo com o princípio da incerteza de Heisen- 
berg. 


E = hu 


( 1 . 2 ) 


A E = hAu 


(1.3) 


A E 


A 


Av = - A - = —1— rí A1Z. (1.4) 

h 2nA t 27r 

onde A E é a incerteza na energia de um nível particular de energia, At é o tempo qne 
um elétron permanecerá no estado de particular energia antes de espontaneamente 
decair para um estado de energia inferior, chamado tempo de vida da transição, e 
Aij é a probabilidade de transição espontânea de Einstein que nos fornece a prob¬ 
abilidade por unidade de tempo de um átomo excitado no nível i decair para um 
nível inferior j com emissão de um fóton. O alargamento natural, em termos do 
comprimento de onda torna-se: 


AA = 


A 2 


— ^ - 


27 TC \At J 


-A, 

27 TC 13 


(1.5) 


1.1.2 Alargamento Doppler 

As partículas que compõem o gás de uma atmosfera estelar, estão em constante 
movimento, e este movimento desloca os comprimentos de onda, visto por um ob¬ 
servador. As linhas espectrais são então alargadas, já que os átomos têm velocidades 
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diferentes e os comprimentos de onda se deslocam pelo efeito Doppler. O alarga¬ 
mento Doppler é o alargamento de linhas espectrais devido ao efeito Doppler em 
que o movimento térmico dos átomos ou moléculas deslocam a frequência aparente 
da radiação emitida. As distintas velocidades das partículas componentes do gás 
provocam pequenos deslocamentos, e o efeito cumulativo é que alarga a linha. Para 
um gás em equilíbrio térmico, alargamento é dependente unicamente da frequência 
da linha, da massa das partículas que irradiam e da temperatura. Portanto o alarga¬ 
mento Doppler pode ser usado como um método para medir a temperatura de ob¬ 
jetos observados. Para velocidades térmicas não relativísticas, o efeito Doppler na 
frequência é dado por: 


v = u Q (l + (1.6) 

onde ué a frequência observada, u 0 é a frequência da fonte de radiação, v é a 
velocidade da fonte de radiação em relação ao observador, ecéa velocidade da luz. 

Como há uma distribuição de velocidades em todas as direções, o efeito cumula¬ 
tivo será o de alargar a linha observada. Para baixas densidades, a distribuição de 
velocidades em qualquer direção é dada pela distribuição de Maxwell. Se P(v)dv é 
a fração de partículas com velocidades entre v e v + dv em qualquer direção, então: 

p{v)dv = \f^m exp (~i w) iv (17) 


onde m é a massa das partículas, T a temperatura e k é a constante de Boltzmann. 

Escrevendo o efeito Doppler em termos da velocidade v — c — 1 j, e substi¬ 
tuindo a velocidade em função do efeito Doppler, tal que, ^ = f- obtemos: 


P{u)du = 



( 


m 


2nkT 


exp 


m 


cD-l 


! \ 


2 kT 


dv 


( 1 . 8 ) 


Podemos reescrever esta expressão como: 


/ mc 2 ( mc 2 (v-v 0 ) 2 

PW)dv = V ímw exp - 


dv 


(1.9) 
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Analisando a equação acima observamos que, 


P{u)du = C exp 


Au 


Au 


Doppler 


( 1 . 10 ) 


onde 


Doppler 

C 


2kT 


m 


(LU) 


é o alargamento provocado na frequência devido ao alargamento Doppler e C é uma 
constante de normalização. Escrevendo a equação acima em termos do comprimento 


de onda e lembrando que |^j = |-^|, obtemos que a espessura da linha espectral 
devido ao alargamento Doppler é dada por: 


A , 2À 2kT 

A\ = —\ . 

c V m 


( 1 . 12 ) 


1.1.3 Alargamento Colisional 

Todos os fenômenos que causam alargamento nas linhas espectrais devido à per¬ 
turbações de partículas vizinhas são chamados de alargamento colisional, ou de 
alargamento por pressão, já que a fonte microscópica de pressão é a colisão de 
partículas com uma parede real ou imaginária. 

O alargamento colisional é devido a colisões de átomos com outros átomos, 
moléculas ou íons. Essas irão perturbar os estados atômicos ou moleculares e levarão 
a um alargamento da linha que é em geral, muitas vezes maior do que a sua largura 
natural. Este aumento na largura da linha é uma função da densidade dos vários 
tipos de perturbadores e é conhecido como alargamento por pressão. Portanto, se 
A Vj(r) e AVi(r) são, respectivamente, as mudanças nas energias dos níveis supe¬ 
riores e inferiores de uma dada transição devido à influência de um perturbador a 
uma distância r de um átomo excitado, a frequência angular da transição sofre a 
seguinte alteração: 


Ato ji 


AVj(r) - AViir) 

h ' 


(1.13) 


Se o potencial é central, isto é, só depende da distância r, é possível representar 
a interação de longo alcance entre o átomo excitado e um perturbador por um 
potencial da forma. 


AVAr) = 


(1.14) 
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onde é uma constante que depende do nível excitado envolvido e também do tipo 
de perturbador. O valor de n está associado ao tipo de perturbador e processo físico 
envolvido da seguinte forma: 

• n — 2, aplica-se ao caso do hidrogênio e íons hidrogenoides no campo elétrico 
externo ou produzido por outros íons ou elétrons. Essa forma de alargamento 
é conhecida como alargamento Stark linear dos níveis de energia, e depende 
linearmente da intensidade do campo elétrico; 

• n = 4, aplica-se ao caso que descreve o alargamento Stark quadrático por 
exemplo do hélio e outros sistemas onde a separação é uma função quadrática 
do campo elétrico dos perturbadores; 

• n — 3, aplica-se ao caso da interação ressonante dipolo-dipolo. Esta interação 
é finita somente quando o átomo excitado interage com um átomo idêntico no 
estado fundamental e quando há uma forte transição de dipolo elétrico permi¬ 
tido, geralmente a linha de ressonância do átomo, que liga os dois níveis. Esse 
alargamento é conhecido como alargamento ressonante ou auto-alargamento, 

• n — 6, aplica-se ao caso de interação dipolo-dipolo atrativa de longo alcance 
que existe sempre entre quaisquer dois átomos. Alargamento desse tipo é 
conhecido como alargamento de van der WaaVs. 

Os dois primeiros efeitos são importantes apenas em gases altamente ionizados 
enquanto que os dois últimos são mais significativos em gases não ionizados. Entre¬ 
tanto, outras contribuições para forças atrativas de longo alcance existem. 

As forças de curto alcance, que são responsáveis pelas ligações químicas ou re¬ 
pulsão interatômica, diferentemente das de longo alcance, não possuem uma ex¬ 
pressão geral, e são de difícil obtenção. No entanto, para o caso de linhas espectrais 
alargadas por colisões com hidrogênio, hélio e neônio, estas forças repulsivas desem¬ 
penham um papel muito significativo. 

A equação (1.14) é obtida a partir do caso de um sistema binário de um átomo 
excitado e um perturbador a uma dada distância e orientação fixa. Mas, em geral, 
o átomo excitado pode interagir com vários perturbadores e é necessário obter uma 
média sobre as distintas orientações e caminhos desses perturbadores. Essa média 
somente pode ser realizada de forma satisfatória em dois casos limites conhecidos 
como aproximação quase-estática e aproximação de impacto ou aproximação de 
deslocamento de fase. O critério que define o intervalo de validade de cada uma 
das aproximações está associado com o tempo de colisão t c . A forma do perfil de 
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linha para uma frequência de separação Acu = cu 0 — cu a partir do centro da linha 
é determinada pelo trem de ondas emitido durante o intervalo de tempo Aí, onde 
Aí ~ A-. Se o átomo emite um trem de ondas até ser perturbado por uma forte 
colisão, o tempo de interesse, Aí, será o tempo médio entre as colisões t c . E uma 
aproximação que pode ser considerada como se o perturbador estivesse ultrapassando 
o átomo emissor em uma trajetória clássica. Para os perturbadores neutros, essas 
trajetórias seriam aproximadamente linhas retas, enquanto que para íons ou elétrons 
seriam hipérboles. Pode-se então considerar que o átomo emissor será perturbado 
somente durante a duração da perturbação, definida por: 

tc = | (1.15) 

onde pé a distância de maior aproximação (parâmetro de impacto) eréa velocidade 
média nesse ponto. 

Quando a duração de uma colisão (í c ) é muito maior do que o tempo entre 
colisões (T c ), a linha alargada é produzida enquanto o perturbador move-se apenas 
uma pequena distância. Portanto, é possível ignorar o movimento dos perturbadores 
completamente e temos a aproximação quase-estática: 

tc > ~ T c . (1.16) 

Essa aproximação é usada para calcular o alargamento Stark produzido por íons 
em um plasma e para alargamento de átomos neutros em pressões > 100 Torr. A 
partir das equações (1.15) e (1.16), observa-se que a aproximação quase estática é 
melhor aplicada em altas densidades, onde T c é menor, e para baixas temperaturas, 
onde v é pequeno. 

No outro extremo, onde a duração da colisão é pequena comparada com o tempo 
entre colisões, a aproximação de impacto pode ser usada: 

íc< Ã^~ Tc ' ( L17 ) 

Na aproximação por impacto o deslocamento de fase produzido por uma colisão 
é calculado e o resultado é uma média sobre todos os parâmetros de impacto. A 
aproximação de impacto é geralmente usada para descrever o alargamento das linhas 
para pressões entre 0 e 100 Torr. Analisando as equações (1.15) e (1.16) observa-se 
que essa aproximação é melhor aplicada para baixas pressões e altas temperaturas, 
onde T c e v são grandes. 

Entretanto, deve ser recordado que o tempo de interesse é determinado pelo 
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deslocamento da frequência a partir da linha central, Ac j. Portanto, é possível que 
partes do perfil de linha sejam calculadas pela aproximação quase estática e outras 
sejam descritas pela aproximação de impacto. 

1.2 Contexto teórico 

Linhas satélites presentes nas asas distantes dos perüs de linhas das séries de Lyman 
e Balmer do hidrogênio atômico têm sido observadas em estrelas e em plasmas 
em laboratório (e.g., Kiclkopf et ah, 2002). Entretanto, a descrição adequada do 
fenômeno encontra dois grandes problemas. 

A primeira dificuldade é do lado experimental, que tem sido o desenvolvimento 
de técnicas para observar as interações do hidrogênio atômico neutro em densidades 
suficientemente altas para efeitos de alargamento das linhas espectrais serem ob¬ 
servados, e para caracterizar as condições quantitativas sob as quais os espectros 
são criados. Por outro lado, do ponto de vista teórico, desenvolver métodos para 
calcular as contribuições a partir de todos os possíveis estados moleculares do H 2 
e H^~ para o espectro de átomos livres na colisão, pois outro absorvedor contínuo 
hidrogenoide é a molécula de hidrogênio H 2 , cuja abundância excede a de hidrogênio 
atômico na atmosfera de estrelas tipo M (T e fj < 2 500 K). Embora essa molécula 
não absorva na região do visível, seus íons H^~ e 14% o fazem. A densidade dos íons 

é proporcional ao produto da densidade de átomos de H neutros pela densidade 
de prótons, portanto ele contribui de forma significativa na opacidade contínua em 
intervalos de temperaturas e pressões dentro das quais existem prótons e átomos 
de ff neutro em quantidades significativas. Esse intervalo de pressão e temper¬ 
atura ocorre precisamente na fotosfera de estrelas tipo A (T e fj ~ 10 000 K), nas 
quais existem aproximadamente as mesmas quantidades de ff neutro e ionizado. 
Consequentemente, aproximadamente 10% da opacidade contínua no visível dessas 
estrelas é devida à contribuição do íon H 2 . A contribuição deste íon para a opaci¬ 
dade contínua foi estudada por Bates (1951), Bates (1952), Bates et ah (1953) e 
Matsushima (1964) entre outros. Matsushima (1964) demonstrou que, para estrelas 
frias, a contribuição do íon é inferior a 10% da contribuição do ff 2 para À > 4000 
À . Doyle (1968) estudou a absorção produzida pela assim chamada quase-molécula 
ff 2 e descobriu que a sua contribuição é significativa na região ultravioleta para À < 
2 500 Ã e estrelas do tipo solar. Enquanto isso, o íon H^" existe principalmente em 
baixas temperaturas, próprias de estrelas M. 

Experimentalmente, linhas satélites características em 1600 Ã e 1405 Ã nas 
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Koester’s T eff =11000K log g=8.25+spSpec—54153-2683-560 spectra 



Figura 1.2: Espectro observado com o SDSS da anã branca SDSS 

J073708.00+411227.45, com sinal-ruído S/R=68 e magnitude g= 15,77, e em ver¬ 
melho um modelo de D. Koester de 11000 K e logg=8,25 que se ajusta, para com¬ 
paração. 

asas de Lyman-a, associadas com transições livre-livre quase-molecular do H 2 e H^, 
foram observadas nos espectros ultravioletas (UV) de algumas estrelas obtidos com 
o International Ultraviolet Explorer (IUE) e com Hubble Space Telescope (HST), 
Greenstein (1980), Koester et al. (1985) e Nelan & Wegner (1985). Os satélites na 
asa vermelha de Lyman-d na região espectral de 905 a 1187 À, foram observados 
com o Far Ultraviolet Espectros copie Explorer (FUSE). 

As estrelas que apresentam linhas satélites na serie de Lyman são as estrelas 
anãs brancas DA ( WDs ), estrelas velhas do ramo horizontal do tipo espectral A, e 
estrelas À Bootis, estrelas de População I com tipo espectral compreendido entre B 
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e F (Allard et al., 2000). As estrelas anãs brancas são o produto final da evolução de 
mais de 95% de todas as estrelas (e.g., Kepler et al., 2007). Elas consumiram todo 
seu combustível nuclear (que depende da massa individual da estrela) e, portanto, 
sua fonte de luminosidade é dominada pela energia interna residual. Seu esfriamento 
via irradiação é lento, devido à pequena área superficial (um raio típico de uma anã 
branca é da ordem de 10 000 km), logo tem um tempo de esfriamento longo, levando 
aproximadamente 10 10 anos para sua temperatura efetiva diminuir de cerca de 10 5 K 
para próximo de 10 3 K. Portanto, as estrelas anãs brancas mais frias são velhas, e 
guardam em sua estrutura o registro histórico da formação estelar e da evolução de 
nossa Galáxia. 

As observações espectroscópicas no ultravioleta (UV) de anãs brancas revelam 
uma forma da linha muito diferente daquela esperada a partir de um simples alarga¬ 
mento Stark, com linhas satélites na asa vermelha de Lyman-/3 próximas de 1078 
e 1060 Ã (Koester et ah, 1996, 1998). A intensidade desses satélites característicos 
e, na verdade, a forma inteira das asas na série Lyman, é muito sensível ao grau de 
ionização da atmosfera estelar e em plasmas em laboratório, pois é a ionização que 
determina a importância relativa do alargamento colisional por íons, ou partículas 
neutras. 

Portanto as linhas satélites podem ser usadas como um diagnóstico da temper¬ 
atura em plasmas onde o equilíbrio de ionização de Saha se mantém. A forma do po¬ 
tencial desempenha um papel dominante na grande diferença das características do 
alargamento quase-molecular (Allard et ah, 1998a). Transições quase-moleculares 
ocorrem quando um fóton é absorvido ou emitido por um átomo enquanto esse 
átomo interage com uma ou mais partículas vizinhas, que podem ser átomos, íons 
ou moléculas (Allard & Kiclkopf, 1982). Em curtas distâncias, a interação entre o 
átomo perturbado e a(s) outra(s) partícula(s) geralmente leva à formação de linhas 
satélites moleculares nas asas das linhas atômicas (Allard et ah, 1998b,a, 1999). 
As estruturas observadas na asa de Lyman-o; em 1 623 Ã e 1405 Ã são devidas 
a absorções das quase-moléculas H 2 e 11%, respectivamente. A intensidade desses 
dois satélites depende fortemente do grau de ionização da atmosfera estelar e conse¬ 
quentemente da temperatura efetiva T eS e gravidade superficial logg (e.g., Holweger 
et ah, 1994). 

A região de formação de linhas espectrais, chamada de atmosfera, em uma anã 
branca, representa menos que 1CP 13 da massa da estrela e é dominada por hidrogênio 
(anãs brancas DA) ou hélio (DB/DO) e algumas vezes com pequenas quantidade de 
elementos pesados (por exemplo, anãs brancas DAZ, DBZ e DOZ) (e.g., Koester, 
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2010). É através da atmofera estelar que a radiação é emitida para o espaço inter- 
estelar, e grande parte das informações são obtidas através de modelos atmosféricos, 
que revelam parâmetros como T e g, log g e composição química. 

A técnica mais bem sucedida para a determinação da temperatura efetiva (T e g) e 
gravidade superficial (log g) de estrelas anãs brancas DA é a comparação dos perfis de 
linha do hidrogênio observados com espectros sintéticos a partir de modelos teóricos 
de atmosfera (e.g., Tremblay & Bergeron, 2009). O sucesso dessa técnica consiste 
no fato de que os perfis de linha teóricos são sensíveis às variações dos parâmetros 
atmosféricos (Wegner & Schulz, 1981) e que as linhas de hidrogênio são as únicas 
presentes nos espectros de estrelas anãs brancas DA, como segue no exemplo da 
figura 1.2. 

A física utilizada para construir os modelos incluem o modelo de gás ideal, e as 
fontes de opacidade radiativa. As abundâncias relativas dos componentes químicos 
podem ser determinadas através da função partição usando a probabilidade ocupa- 
cional de Hummer & Mihalas (1988), que resolvem a população de cada nível do gás 
com a técnica de minimização de energia livre. 

É necessário conhecer as transições atômicas para explorar ao máximo a técnica 
espectroscópica. Bergeron et al. (1992) foram um dos primeiros a aplicar esta aprox¬ 
imação a um grande número de anãs brancas; eles analisaram as linhas de Balmer. 
A técnica foi mais tarde aplicada às linhas UV da série de Lyman (Kepler & Nelan, 
1993, Bergeron et ah, 1995, Vennes et ah, 2005). A análise espectroscópica é limi¬ 
tada principalmente pela qualidade das observações, tais como a razão sinal/ruído, 
calibração de fluxo, etc. O valor absoluto dos parâmetros atmosféricos depende dire¬ 
tamente do nível de sofisticação da física incluída nos cálculos do modelo de atmos¬ 
fera, da consistência no cálculo do efeito Stark e da teoria de convecção (Tremblay 
& Bergeron, 2009). 

Em 2002 Bergeron & Leggett (2002) publicaram melhorias nos modelos, levando 
em conta o efeito nas opacidades óticas quase-moleculares dos momentos de dipolo 
induzidos por interações entre as partículas. Essa interação provoca um rearranjo 
temporário das cargas eletrônicas, com formação de uma molécula temporária, pro¬ 
duzindo assim absorções ou emissões além daquelas dos átomos isolados. 

Atualmente existe um grande número de espectros óticos de anãs brancas obtidas 
como parte do Sloan Digital Sky Survey (SDSS) (Eisenstein et al., 2006, Kleinman, 
2010), permitindo assim o estudo da distribuição de massa das anãs brancas (Kepler 
et al., 2007). Kleinman et al. (2004) e Liebert et al. (2005) já tinham estimado as 
massas a partir de análise espectroscópica. Bergeron et al. (1991) e Koester (1991) 
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discutem o aparente aumento na massa média de estrelas anãs brancas DA com 
Tcflf <12 000 K. Kepler et al. (2006) compararam os espectros do SDSS com espec¬ 
tros de alto sinal/ruído (S/R) obtidos com o telescópio Gcmini 8 m para quatro 
anãs brancas com T e fj ~ 12 000 K e determinaram que as massas derivadas a partir 
da análise dos espectros do SDSS, com S/R relativamente baixos, estão sistemati¬ 
camente superestimadas por AM ~ O.13M 0 . Devido à correlação entre T e s e log g 
(um pequeno aumento em T eS pode ser compensado por uma pequena diminuição 
em logg), tais discrepâncias estão concentradas unicamente nas regiões em torno 
do máximo das linhas de Balmer, 14000 K > T e g > 10 000 K. Para estudar a 
tendência de aumento da massa média aparente das anãs brancas DAs, Kepler et al. 
(2010) analisaram 1505 delas com S/N > 20 e T e fj > 12 000 K e determinaram 
uma massa média das DAs de M/M & = 0.604 ± 0.003. Eles observaram que a 
distribuição da amostra é similar à publicada por Liebert et al. (2005) obtida do 
Palomar Green Survey com 348 estrelas. Falcon et al. (2010) usaram uma média 
de conjunto de avermclhamento gravitacional de 449 DAs, método independente 
dos perfis de linha pois usa somente a posição do centro da linha, e estimou uma 
massa média de 0.640 ± O.O14M 0 , e que não apresentava mudanças significativas 
para T e fj < 12 000 K. O valor de massa média obtido por Falcon et al. (2010) é 
similar ao obtido por Tremblay & Bergeron (2009) para as estrelas quentes quando 
incluem o formalismo de Hummer & Mihalas (1988) para o cálculo do campo elétrico 
médio causador do efeito Stark. 

Koester et al. (2009), Koester (2010), Tremblay & Bergeron (2009) e Gianninas 
et al. (2010) estudaram as possíveis causas para o aumento aparente na massa deter¬ 
minada espectroscopicamente para T e g < 12 000 K, não observado nas determinações 
fotométricas. Eles concluíram que esse fenômeno não é causado por contaminação 
de He, e Tremblay et al. (2011) propôs que o uso de modelos unidimensionais nos 
cálculos da convecção da atmosfera, comparados a um modelo mais realístico tridi¬ 
mensional, poderia ser a causa. Kepler et al. (2010) descobriram campos magnéticos 
acima de 1 MG em cerca de 5% das anãs brancas, e que o campo magnético tende a 
ser maior quando a estrela desenvolve uma zona de convecção atmosférica profunda, 
em T e g < 12 000 K e propõe que uma causa do possível alargamento de linha obser¬ 
vado, interpretado como aumento de gravidade superficial, possa ser alargamento 
magnético (efeito Zeeman) para campos abaixo do limite de detecção observacional. 
Há um enorme esforço para eliminar as inconsistências nos cálculos do modelo, e es¬ 
sas tentativas mostram a necessidade do contínuo estudo da interação do hidrogênio 
perturbado por prótons, elétrons, átomos e moléculas em alta pressão. 
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As estrelas anãs brancas possuem uma alta densidade em sua atmosfera, o que 
torna o alargamento colisional (alargamento por pressão) dominante em relação aos 
demais. As colisões ou interações entre partículas e um átomo perturbado induzem 
mudanças nos níveis de energia dos átomos perturbados e correspondentes alterações 
nas frequências das transições atômicas. Essas colisões entre átomos perturbados 
e partículas perturbadoras são responsáveis pela formação temporária de quase- 
moléculas que levam ao surgimento de possíveis satélites característicos nas asas do 
perfil de linha, ou seja, absorções ou emissões além das esperadas. A figura 1.3 
mostra o espectro da anã branca G226-29, observada com o Hubble Space Tclescope 
(HST) por 6 exposições, cada tempo de 3 horas, entre setembro de 1994 e dezembro 
de 1995 por Keplcr et ah (2000). O espectro da G226-29 mostra a existência e forma 
dos satélites. 


G226-29 Average HST Spectra 
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Figura 1.3: Espectro observado com o Hubble Space Telescope (HST) da anã branca 
G226-29, de Kepler et ah (2000), mostrando a asa vermelha de Lo:, centrada em 1 215 
Ã e as linhas satélites em 1405 Ã do e 1 600 Ã, do H 2 

Portanto, o perfil de linha de estrelas anãs brancas, devido a interações entre o 
hidrogênio atômico e prótons, é distinto daquele calculado levando em consideração 
apenas partículas isoladas da atmosfera. Há a formação da quase-molécula II 2 e sua 
dissociação. E quando a diferença AH( R ) entre um potencial de um nível superior 
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e um potencial de um nível inferior apresenta um extremo, haverá a formação de 
um satélite centrado periodicamente em Acu = kAV ext (R), com (k = 1,2,3,...) 
(Anderson, 1956, Allard & Kielkopf, 1982). A variação da frequência Acu é medida 
a partir do centro da linha espectral não perturbada e o satélite característico, e 
escrita na mesma unidade da diferença de potencial, por conveniência. A posição de 
A Vext(R) e a dependência de AV(R) com a separação internuclear R determinam a 
amplitude e a forma dos satélites, respectivamente. 
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1.3 Visão geral da tese 

Nessa tese apresentarei o estudo da interação do átomo de hidrogênio com prótons, 
que temporariamente cria a molécula de H^~, com abundâncias mostrada na figura 
1.4, para densidade inicial de 1,0 x 10 17 moléculas de H 2 por cm -3 e pressão con¬ 
stante. 



Figura 1.4: Abundância relativa das quantidades de n(H 2 )/n(H 2 ) em preto e 
n^H^/n^H) em vermelho, com densidade inicial de 1,0 x 10 1 ' H 2 cm -3 . 

No capitulo 2 serão apresentados os cálculos do potencial do H^~ até o nível 
n — 10 (número quântico principal da molécula), bem como a análise de todas as 
possíveis transições que podem ocorrer. Apesar da molécula de ser muito bem 
conhecida, seus potenciais não possuem publicações até o nível n = 10, nem o grau 
de resolução que será apresentado. Esses novos dados teóricos devem ser inseridos 
no modelo de perfis de linha para altas densidades, em adição ao alargamento Stark. 

No capítulo 3, apresentarei nosso estudo dos momentos de dipolo elétrico para 
todas as possíveis transições. Da mesma forma que os potenciais, na literatura 
encontramos os valores somente para Lya, Ly/3, Lyy e Ha, sendo que os momentos 
de dipolo elétrico para Lyy foram usados por Allard et al. (2004), mas não foram 
publicados até o momento. 

O capítulo 4 é destinado aos cálculos dos perfis de linhas das transições das 
séries de Lyman e Balmer, bem como a análise das linhas satélites e quais transições 
produzem tais linhas. 








Introdução 


17 


Os capítulos 3 e 4 foram publicados por Santos & Kepler (2012) no Monthly 
Notices of the Royal Astronomical Society, Volume 423, Issue 1, pp. 68-79, em 
junho de 2012, DOI: 10.1111/j. 1365-2966.2012.20631.X 

Finalmente, no capítulo 5 serão apresentadas as discussões e conclusões referentes 
à essa tese. 




Capítulo 2 

Potenciais da molécula de 
Hidrogênio H(}~ 


A estrutura do íon molecular foi estudada por Teller (1930), Hylleraas (1931), 
Jaffé (1934), Baber & Hassé (1935), Bates & Poots (1953) e Madsen & Peek (1971). 
Entretanto, tabelas precisas da energia eletrônica em função da distância internu- 
clear para estados excitados n > 3 não foram até agora publicadas. 

Como originalmente demostrado por Burrau (1927), a equação de Schrõdinger 
para o H:)" pode ser separada e a energia eletrônica calculada. 

2.1 A aproximação de Born-Oppenheimer 

As funções de onda e energias de uma molécula são obtidas a partir da equação de 
Schrõdinger: 


HV(R,r) = EV(R,r) (2.1) 

onde R e f representam as coordenadas nucleares e eletrônicas respectivamente. 

Considerando uma molécula diatômica composta por núcleos e elétrons com 
massas pontuais e desprezando as interações de spin-órbita e outras correções rela- 
tivísticas, o Hamiltoniano molecular é dado por (Levine, 2004): 


77 V \ ' 1 Y7 2 \ ' 7 r 7 , \ ' \ ' ZaZpe 2 


2 z ' m 

a 

•v 

f N (R) 


O-^a/3 


A(rO 


47re n A, 

a p>a 

v 

V NN (R) 


\ ^ \ ^ Z a e 2 
Z-J 47T e 0 r ic 


EE 


j *>i 


47re 0 r, 




V eN (r,Ê) 


V e e{r) 

( 2 . 2 ) 
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onde a e (3 se referem aos núcleos e i e j aos elétrons. 

Os cinco termos do Hamiltoniano (2.2) representam os operadores: 

Tn(R) energia cinética dos núcleos; 

T e (r) energia cinética dos elétrons; 

Vjvjv(-R) repulsão entre os núcleos, sendo R a g a distância entre os núcleos a e (3 
de números atômicos Z a e Zp | 

V e N r(r, R) interação entre os elétrons e os núcleos, onde r ia é a distância entre o 
elétron i e o núcleo a, 

V ee (f) interação entre os elétrons, onde é a distância entre os elétrons i e j. 

O termo VAv('C R) infelizmente impede a separação do H em componentes nu¬ 
clear e eletrônica, que permitiria escrever a função de onda molecular como um pro¬ 
duto de termos nuclear e eletrônico, \E f ( R,r ) = x(A)0(r). Torna-se então necessário 
utilizar algum tipo de aproximação. Introduzimos a aproximação de Born-Oppenheimer. 
O termo V eN (r,R) não pode ser desprezado, mas pode-se parametrizar R de modo 
que a função de onda total possa ser reescrita como x(f?)0(r; R). A aproximação de 
Born-Oppenheimer consiste no fato de que os núcleos são muito mais massivos que 
os elétrons, permitindo considerar que os núcleos estão fixos enquanto os elétrons 
realizam o movimento. Fixando R , pode-se resolver a função de onda eletrônica 
(f){r\ R) com dependência apenas parametricamente em R. Com isto, obtém-se a 
curva de energia potencial ao longo do movimento parametrizado do núcleo. Con¬ 
sequentemente, a função de onda pode ser reescrita: 

V(R, r) = x{R)(t>(?\ R). (2.3) 

e o Hamiltoniano, em uma forma simplificada e em unidades atômicas (Apêndice 
A): 

è = £ v ? + nr-,Ã) (2.4) 

a i 

onde V (r, R) é a energia potencial total do sistema. 

Substituindo 2.3 e 2.4 em 2.1 obtém-se: 

\r\ E — v « - 5 E + V ( F ’ «)1 x(R)m r) = E X (R)m r) ( 2 . 5 ) 

a. i 
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ou seja, 

“E —+ Í-^E V 2 i + v (f,R)}x(R)<t>(?-,R) = Ex(R)(/>(r; R), 

Z 777-q- L J l Z 

a i 

(2.6) 

sendo 

V„ [x W(f; £)] = ximlttr; R) + V0(r; R) ■ V X (R) + 0(r; R)V 2 aX {R) 

onde 'S7(j){f]R ) • W X (R) representa o acoplamento elétron-núcleo. 

De acordo com a aproximação de Born-Oppenheimer, Vd(r: i?) « 0 e V^0(r; i?) ~ 
0 , pois a nuvem eletrônica reorganiza-se quase instantaneamente para qualquer 
mudança na posição nuclear, significando que uma pequena perturbação na con¬ 
figuração nuclear é desprezível, motivo pelo qual essa aproximação também é con¬ 
hecida como aproximação adiabática. Resta então: 

V* [x W(r; R)] = 0(r; R)V 2 aX (R) 

que é aplicada em 2.6: 

R) E —' Vlx(R) + H E V * + V (e £)1 x(R)<t>(r; R) = E X (R)<t>(f-, R), 

Z 777-0; L Z 

a i 

(2.7) 

Devido, ainda, à aproximação de Born-Oppenheimer, o primeiro termo, corre¬ 
spondente ao operador energia cinética nuclear da equação 2.7, é desprezado, e o 
termo que descreve a repulsão nuclear [Vjvjv(-R)] torna-se uma constante e pode, 
sem perda de generalidade, ser adicionado ao hamiltoniano, sem alteração de suas 
autofunções. Portanto, pode-se dividir ambos membros da equação 2.7 por x(R) 
resultando em: 

E V " + V(f, R)] 0(r; R) = E e (j)(f; R), (2.8) 

1 

onde E e é a energia eletrônica do sistema para uma dada configuração nuclear. 
Consequentemente a equação de Schrõdinger eletrônica pode ser escrita como: 

H e (f>(f- R) = E e (j>(r :; R) (2.9) 

A obtenção da função de onda nuclear depende diretamente da resolução da 
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equação 2.8, ou seja: 


■E(r;Ã)£— V ^(fl) + 

2 ' rn a 


ij2^+nr,R)+j:j: ZaZf 


a (3>a 


R, 


■a (3 


x(R)<P(r;R) 


= Ex{R)<t>{r .; R) 

( 2 . 10 ) 


Simplificando a 2.10 e usando a 2.8, obtém-se: 


1 

2 


E — V «X(^) + [^e + EE 
L ~ 

a p>a 


Z a Zp~ 

Raft - 


x(R) = E X (R) 


( 2 . 11 ) 


O termo entre colchetes da 2.11 é o potencial efetivo do sistema e a grandeza 
E é a energia total, que envolve as energias eletrônica, vibracional, rotacional e 
translacional da molécula. 


2.2 O íon molecular do Hidrogênio H 2 

A estrutura do H 2 é formada por um elétron e dois prótons, conforme a figura 2.1. 
O problema do íon molecular do ffidrogênio é de grande interesse, tanto para físicos 
quanto para astrofísicos, pois o e outras moléculas diatômicas com um elétron 
são as únicas entre moléculas em que a equação de Schrodinger para seu movimento 
eletrônico pode ser separável em três equações diferencias ordinárias (Madsen & 
Peek, 1971), devido ao fato que os núcleos são muito mais massivos e lentos em 
relação ao movimento eletrônico, ou seja, na aproximação de Born-Oppenheimer. 



Figura 2.1: Estrutura do íon molecular do ffidrogênio H^ [Adaptado de Levine 
(2004)] 
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Portanto, pode-se escrever a partir da equação 2.8, a equação de Schrõdinger da 
parte eletrônica, que possui a seguinte forma em unidades atômicas: 



Na figura 2.1 a origem de coordenadas está no centro do segmento de reta que 
une as cargas positivas e essa linha coincide com o eixo z em coordenadas carte¬ 
sianas. Devido à simetria do sistema apresentado nessa figura, tende-se a resolver, 
em primeira instância, a equação de Schrõdinger eletrônica utilizando as coordenadas 
polares esféricas (r, 9, 0) do elétron, ou seja, escrevendo r a e r b nessas coordenadas. 
Entretanto, utilizando as coordenadas polares esféricas, a equação de Schrõdinger 
não pode ser separada. 

Por sua vez, se forem utilizadas coordenadas esferoidais prolatas, também con¬ 
hecidas como coordenadas elípticas confocais (À, /i, 0), conforme figura 2.2, é possível 
fazer a separação de variáveis e obter soluções, como demonstrado por Burrau (1927). 
As coordenadas elípticas são usadas em problemas envolvendo dois centros, a e b, 
separados por uma distância fixa R, Se r a é a distância do núcleo a até a posição 
do elétron, e r b é a distância do núcleo b até o mesmo ponto, então as coordenadas 
elípticas são definidas como: 


_ T a + n 

R 


(2.13) 


ú 


r a ~ r b 

R 


(2.14) 


0 — tan 1 j (2-15) 

À define clipsóides de revolução com focos a e b e n define parabolóides de revolução 
em torno do eixo z, 0 é o ângulo de rotação do elétron em torno do eixo internuclear 
z, o mesmo usado em coordenadas polares esféricas. 

Os intervalos dessas coordenadas são: 


1 < À < 00, 

-1 < n < 1 , 

0 < 0 < 27r. 

As coordenadas elípticas da figura 2.2 estão relacionadas com as coordenadas 
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Figura 2.2: Coordenadas esferoidais [Adaptado de Li et al. (2001)] 


cartesianas pelas seguintes transformações: 


X = — COS 0a/(/ 2 - 1)(1 - A 2 ) 

(2.16) 

y = ^sen4 V(A 2 - 1)(1 - ^ 2 ) 

(2.17) 

RX/i 

2 

(2.18) 

e 0 elemento de volume é dado por: 


r>3 

dr = —(X 2 — n 2 )dXdjjíd(j). 

(2.19) 


A equação de Schrõdinger 2.12 pode ser escrita agora em função do novo sistema 
de coordenadas. Para tanto, é preciso obter r a , r b e o V 2 em coordenadas elípticas. 
Os valores de r a e r b são obtidos diretamente das equações 2.13 e 2.14, fazendo a 
soma e subtração de ambas: 
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r a = ^(A + /i) 


( 2 . 20 ) 


n = 


( 2 . 2 !) 


A expressão do Laplaciano em coordenadas elípticas é obtida a partir do Lapla- 
ciano em coordenadas curvilíneas ortogonais, pois o sistema de coordenadas esfer- 
oidais prolatas é um sistema curvilíneo ortogonal. O Laplaciano em coordenadas 
curvilíneas é dado por (Kaplan, 1953): 


V 2 = 


d_ í h^hçj, d \ d_ í h^hx d 

h\hfihff) [< 9 A \ h x dXj dfi \ h tJ d\i 


d_ (hxhft d 

r\ 1 


d(f) V K d(j)J \ 


( 2 . 22 ) 


E seus coeficientes métricos são dados por: 


fdx \ 2 (dy \ 2 fdz \ 2 R /A 2 — // 


\dXJ 


hx — \ \ TTT I + ^ ^ 


d\ 


d\ ) 


A 2 


(2.23) 


, fdx \ 2 fdy \ 2 (dzy 

='/ d ) + d) + d> 


R I A 2 -/i 2 

2 11^7 


(2.24) 


Substituindo 2.23, 2.24 e 2.25 em 2.22 obtém-se o Laplaciano em coordenadas 
elípticas: 


V 2 = 


R 2 (A 2 - y 2 ) {dX 


d 


d 


C ~ 4 à 


d_ 

ô/i 




+ 


(A 2 - d) 


d 2 


(A 2 — 1) (1 — y 2 ) d(j) 2 
(2.26) 

Com as equações 2.20, 2.21 e 2.26 reescreve-se a equação de Schrõdinger eletrônica 
em coordenadas elípticas da seguinte forma: 


R 2 (X 2 - y 2 ) 


d 

dX 


»■-< 


d_ 
^ dy 
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(A 2 - y 2 ) 


d 2 


(A 2 — 1)(1 — /i 2 ) def ) 2 




M (v 1 . T 

-= E -vb 

R{X 2 -y 2 ) V R 


(2.27) 
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r 1 -1 

Multiplicando a equação 2.27 por — R2 ^2 e usando a notação £ = E — 

para a energia eletrônica, temos 


2 . í ÍA 2_ ,P 1 , r, _ 2l _a 1 (a 2 - p) g \* 

0A L SAj 0M L ayjJ (A 2 - 1)(1 - p d<t> 2 J 


eR 2 (X 2 — /i 2 ) 

+2AÍM =---— V 

2 


(2.28) 


Por conveniência, define-se um parâmetro (real) p via p 2 = £ 2 R e reescreve-se a 

equação 2.28: 

JA U _ n Al , A [n _ ^ A1 + A ~ A d 2 \ * , 

\<9A |_ ^ + (A 2 -l)(l-/x 2 )d0 2 J + 


+2\m = p 2 (X 2 - p 2 )^ 


(2.29) 


passando todos os termos para o lado esquerdo, 


A [ fA 2_ n A'1 + A |m_ u 2 ) Al + JAiíLAI * , 

<9A 1 ôaJ + dp |_ ^ ’dp\ + (X 2 p 2 )d(j) 2 j + 


+ [2XR - p 2 (X 2 - p 2 )] 'P = 0 


(2.30) 


A equação (2.30) é completamente separável, resultado já esperado pelo uso das 
coordenadas elípticas, e isso permite que a função de onda seja escrita como 


tt(A,/^) = A(A)M(//)*(0) (2.31) 

Fazendo a substituição de ^(A ,p,(f>) = A(A)M(/i)$(0), na equação 2.30, obtém- 
se a equação separada em três equações diferenciais ordinárias (Bates et al., 1953), 

’ d 2 

— + m 2 J S = 0 (2.32) 

f [(! - pm + \-A+pY - -^J M = 0 (2.33) 

op l op J |_ 1 — p- 

A (A 2 -d| + A + 2RX-p 2 X 2 -^- í A = 0 (2.34) 
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onde A e m 2 são as constantes de separação. As constantes m, A e p são determi¬ 
nadas pelas condições de que a (A, /q 4>) seja única, contínua e finita sobre toda a 
região do espaço definido por 1 < À < 00 , — l<p<leO< 0 < 2ti. 

A primeira equação, equação (2.32), refere-se à parte azimutal do problema en¬ 
quanto que as equações (2.33) e (2.34) referem-se às partes angular e radial, respec¬ 
tivamente. 

A equação azimutal (2.32) pode ser resolvida diretamente com uma solução da 
forma: 

$( 0 ) = e (± im ú) ( 2 . 35 ) 

para que a função de onda tenha solução única, onde m, a constante de separação, 
deve assumir valores inteiros (m = 0,1,2,3,...) e está associada à projeção do 1110 - 
mentum angular orbital do elétron sobre o eixo nuclear. A designação do orbital 
está associada ao valor de m, ou seja, a componente o corresponde a m — 0 , e a 
componente 7 T a m — 1 , e assim por diante. 

2.2.1 Função angular 

A solução da parte angular 2.33 pode ser obtida a partir do método sugerido por 
Hylleraas (1931), que consiste em representar a solução de 2.33 como uma expansão 
em polinómios associados de Legendre na forma: 

OO 

M(íO = X>J7» (2-36) 

l=m 


onde P; m (p) é o Z-ésimo polinómio de Legendre, solução da equação diferencial de 
Legendre, 



(1 - » 2 ) 


8PT M 

drj 


+ 


1(1 + 1 ) 


9 

m 

1 — rj 2 


PT(v) = 0 


(2.37) 


que possui a mesma forma da equação 2.33. Substituindo 2.36 em 2.33, a parte 
angular toma a forma: 


l=m 



(1-V) 


dirivY 

dn 


+ 


2 ,.2 


-A + p n 


2 

m 

1 — /i 2 



0. (2.38) 
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Fazendo agora a substituição de 2.36 na equação diferencial de Legendre 2.37, obte¬ 
mos: 


OO 

2> 

l=m 


d_ 

<9/j 


(1 


2 ,9<rw 


ô/i 


-/(/ + 1 ) - 


m 




P/» = 0 


A relação entre a parte angular e a equação de Legendre é direta, pois 
a primeira parte do lado esquerdo da equação angular 2.38 temos: 


(2.39) 

isolando 


EU£ 


l=m 


(9/i 


( 1 -^) 




(9/í 


OO 

= E« 

l=m 


+A - p 2 /i 2 + 


m 


l-/i 2 j 


PTb 0 ( 2 . 40 ) 


Substituindo na equação de Legendre, encontramos: 


E 

l=m 


Cl 


m 


+A — p fi + --- — 1(1 + 1) — 

1 — /i- 2 


m 


1 — /i 2 


pnti = 0 (2.4i) 


ou seja, 


OO 

E c ' l+ A - pV -'(' +!)] ÍTM = 0 

l=m 

Usando a relação de recorrência dos polinómios de Legendre P ; m (/i), 

(/ + mJ-P^Qi) + (f ~ m + l)P/+i(/x) 

2 / + 1 

duas vezes para eliminar os termos /i e /i 2 obtemos 


/iP/» 


(2.42) 


(2.43) 


/i 2 P”>) = 


(/ + m — 1 )(/ + m) 


L ( 2 /-l)( 2 Z + l) J 

(Z — m + l)(l + m + 1 ) | (l — m)(l + m) 

( 2 / + 1)(2Z + 3) + ( 2 / + 1 ) (2Z — 1) 

(l — m + 2)(l — m + 1 ) 
+ (2Z + l)(2/ + 3) 

Substituindo o termo 2.44 na equação 2.42, temos 


+ 

pr (/i) + 


]T Q [A-/(/ + l)]P™(/i)-p 2 

l=m 


(l + m — 1 )(/ + m) 

(21 — 1 )( 2 / + 1 ) 


%(/*) + 


(2.44) 
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-P 


(l — m + l)(l + m + 1) (/ — m)(l A m) 


(2/ + l)(2/ + 3) (2/ + 1)(2Z — 1) 

2 (l — m A 2) (7 — m A 1) 


pr(p) + 


-p 


(21 + 1)(2Z + 3) 

Agrupando os termos em função dos polinómios de Legendre, 


PllM = 0 (2.45) 


£ 

l=m 


-p 


(l + m — 1 )(/ + m) 
( 2 /-l)( 2 / + l) 


A-l(l + l)-p 


2 1 (l — m A l)(l A m + 1) i (/ — m){l + m) 


+ 


(21 + 1)( 2 Z + 3) (21 + 1 )( 2 Z — 1)/ J 

2 (Z — m + 2)(l — m + 1 ) 


-P 


(21 + 1) (21 + 3) 


ciP^M + 

ciPr(p) + 
ciPr +2 (p) = 0 


(2.46) 


Baixando os índices dos polinomios de Legendre de (2.46) 


+Q % A — 1(1 + 1 ) — p 2 


(l — m — 1)(/ — m) 
Cl ~ 2 ^ [ (2/ — 1)(2Z — 3) J 

(/ — m + 1 )(/ + m + 1 ) | (l — m)(l + m) 
(2/ + l)(2/ + 3) + (21 + 1 )(2Z — 1) J 

2 (/ 4“ 777. + 2)(Z + 171 Al) 

(21 + 5) (2/ + 3) 


+C7+2 < -p 


+ 

+ 

= 0 (2.47) 


«iQ-2 + biCi + d z c /+ 2 = 0 


(2.48) 


A equação (2.48) é a relação de recorrência para os coeficientes q em função dos 
termos cq, 6 ; e df. 


ai = -p 


(l — m — 1 )(/ — m) 
(2/ — 1)(2Z — 3) 


(2.49) 


bi = A — 1(1 + 1 ) — p 2 


(/ — m + 1 )(/ + m + 1 ) (/ — m)(/ + m) 


(21 + 1) (2Z + 3) 


(21 + 1 ) (2Z — 1) 


(2.50) 


di = -p 


2 (/ A fxi A 2)(/ A Txi Al) 


(2.51) 


(2/ + 5) (2/ + 3) 

Os coeficientes da relação (2.48), são obtidos a partir da relação em fração 
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contínua da forma: 


Q 


aiCi-2 + hei + dici+ 2 = 0 
bici + diCi+2 = —aici- 2 
6>Z + dl ° l+2 ') = —aiCi-2 


Cl 

Cl = 


—aiCi-2 

(k + 


(2.52) 


A relação em fração contínua (2.52) possui infinitamente muitas raízes que podem 
estar correlacionadas com as autofunções, possuindo 0, 1, 2, 3, ... nodos no intervalo 
— 1 < H < 1. A obtenção dos coeficientes q, via fração contínua (2.52), ocorre por 
sucessivas aproximações, e no caso limite de l quando l —> 00 , o termo —* 0, e 
Co = 1. Portanto a expansão (2.36) converge no intervalo —1 < p < 1 unicamente 
se m, A e p satisfazem a relação em fração contínua (2.52) e consequentemente os 
coeficientes q devem ser nulos ou mínimos para l —» 00 . 

A relação de recorrência (2.48), da solução da parte angular, é um conjunto in¬ 
finito de equações lineares homogêneas e para ter solução distinta da trivial, o deter¬ 
minante da matriz formada pelos coeficientes deve ser nulo. Observa-se na relação de 
recorrência que somente os termos em / que diferem por 2 estão relacionados, então, 
somente os termos pares ou impares são relevantes para determinado estado. Isso é 
devido ao número quântico l, reponsável pela paridade do orbital. Com isso, tere¬ 
mos duas matrizes, uma para l par e outra para / ímpar e ambos os determinantes 
dessas matrizes precisam ser zero para que a solução seja distinta da trivial. 

Para / par: 


60 djQ 0 0 0 0 

a 2 b 2 d 2 0 0 0 

0 CL4 64 d 4 0 0 

0 0 06 &6 ^6 0 


1 

Co 


" 0 " 


C 2 


0 


c 4 


0 


C 6 


0 


C 8 


0 






(2.53) 
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Para l ímpar: 


61 d \ 0 0 0 0 

a 3 ò 3 d 3 0 0 0 

0 05 Ò5 CÍ5 0 0 

0 0 Q'7 67 CÍ7 0 


1 

Cl 


" 0 " 


C 3 


0 


C 5 


0 


C 7 


0 


C 9 


0 






(2.54) 


Portanto, podemos reescrever a equação (2.36), a solução da equação (2.33), para 
levar em conta de forma explícita nos cálculos o fator de paridade gerado por /. 


M(,,) = 5>,/C + >), (2.55) 


onde o somatório se estende somente sobre valores pares ou ímpares de s. 

O cálculo do determinante da matriz dos coeficientes (par ou ímpar) da relação 
de recorrência (2.48) gera uma função fi(A,e,R), que depende de A, e e R, que 
sozinha não possui solução nem analítica nem numérica, necessitando da solução da 
parte radial para formação de um sistema. 


2.2.2 Função radial 

Para resolver a equação radial (2.34) considera-se a solução na forma proposta por 
Jaffé (1934), isto é: 

A(A) = (A 2 - 1)t(A + iye^y(0 (2-56) 


com 


R 

cr = - m — 

V 


1 


(2.57) 


(A — 1) 
(A + 1) 


(2.58) 


Essa é a modificação da expansão original proposta por Hyllcraas (1931), em 
termos dos polinómios associados de Laguerre de argumento (À — 1), que é matem¬ 
aticamente mais satisfatória. Substituindo a solução (2.56) na equação 2.34, (e.g., 
Bates & Poots, 1953): 


C(! - C) 2 ^ + [~(rn + 2a - 1)( 2 + 2(cr - 2 p 


l)C + m + l]^ + 
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+ [cr(cr + m)( + er(l + 2 p) + m(m + 1 + cr) + A — p 2 ]y = 0 (2.59) 

Expandindo y(Ç) em série de potência, 

OO 

3/(0 = X>C* (2.60) 

t =0 

e substituindo em 2.59, obtém-se a relação de recorrência 

<x t gt +1 - Pt9t + Itdt-i = 0 (2.61) 

com 

a t — (t + l)(í + m + 1) (2.62) 

Pt = 2 1 2 + (4 p — 2a)t — A + p 2 — 2 pa — (m + 1 )(m + cr) (2.63) 

7 1 — (t — 1 — cr) (t — 1 — cr — m) (2.64) 


Os coeficientes da série YltL 0 9tC s ã° obtidos a partir de uma nova fração contínua 
da forma: 


Oí t 9t +1 - Pt9t + lt9t -1 = 0 
a t9t+\ — Pt9t — ~lt9t -1 
9t ( -Pt + —— ) = ~lít9t-i 
—lt9t -1 

9t — / \ 

(-A + ^ 


(2.65) 


A fração contínua (2.65) possui uma infinidade enumerável de raízes correspon¬ 
dendo as autofunções possuindo 0, 1, 2 , 3, ... nodos no intervalo 1 < À < 00 . A 
partir da expressão (2.65) observa-se que <?_i = 0 e quando t tende ao valor extremo 
da expansão (t —> 00 ) e portanto —>• 0 . 

A expressão 7 ^- 1 — (3tgt+®tgt+i = 0 da solução da parte radial também constitui 
um conjunto infinito de equações lineares homogêneas, que para ter solução distinta 
da trivial, requer que o determinante da matriz formada pelos coeficientes seja nulo. 
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— (3 0 a 0 0 0 0 

71 —{ 3 \ «i 0 0 

0 72 -P2 a 2 0 

0 0 73 - p 3 a 3 

O cálculo do determinante da matriz dos coeficientes da relação de recorrência 
(2.61) gera a função f 2 (A,e, R), que depende de A, e e R, que juntamente com a 
função fi(A,e,R), forma um sistema de equações não lineares que somente possui 
solução numérica. Dessa forma encontra-se os valores dos parâmetros A, e e R. 

2.2.3 Solução geral da parte eletrônica 

A solução geral é obtida a partir das soluções da parte azimutal, angular e radial da 
equação diferencial (2.30) e possui a forma: 


0 

0 

0 

0 


1 

9o 


" 0 " 


9i 


0 


92 


0 


93 


0 


94 


0 






( 2 . 66 ) 


' 00 / \ _ 1 \ t 

*(A, n, <t>) = Ne^HX 2 - l)í m (A + 1 ré-r» T T 9, (2.67) 

8 t =0 \ ' / 

onde N é a constante de normalização e é obtida a partir de: 

/»00 /*1 /*27T 

/ / / l^(A,/i,(/>)| 2 du = 1 (2.68) 

J i J- iJo 

onde dv é o elemento de volume. 

Analisando a solução geral, bem como as equações que a geraram, observa-se 
a dependência direta de £ com m 2 , o que origina uma dupla degenerescência para 
valores com m ^ 0. Desta forma existem dois estados possíveis com o mesmo valor 
de m 2 , com números quânticos dados por m e — m e que podem ser duplicados se 
for levado em conta o spin eletrônico e — Há também a dependência com /, 
número quântico responsável pela paridade. Portando cada estado possui muitos 
autovalores que o distinguem dos demais, o que leva a um conjunto de terminolo¬ 
gia dos mesmos. A terminologia dos orbitais moleculares (OM) de uma molécula 
diatômica homonuclear pode ser dada conforme a notação de átomos separados, ou 
notação de único átomo, ou ainda notação de simetria (Levine, 2004). A notação 
de átomos separados leva em consideração o estado do átomo de hidrogênio obtido 
através da dissociação da molécula para grandes valores de R. A notação de único 
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Tabela 2.1: Terminologia dos orbitais moleculares (estados eletrônicos) de uma 
molécula diatômica homonuclear_ 


Átomos Separados 

✓ r 

Unico Atomo 

Simetria 

o g ls 

lsa g 

la g 


2 P< 

1 <J U 

a g 2s 

2sa g 

2a g 


3 P< 

2 cr u 

n u 2p 

2pvr.u 

iTTu 

a g 2p 

3 SUg 

3 <Tg 


Tabela 2.2: Representação de l. 


Ipar 

letra 

limpar 

letra 

0 

s 

1 

p 

2 

d 

3 

f 

4 

9 

5 

h 


átomo, também chamado de átomo unido, leva em consideração o átomo formado 
quando a distância internuclear tende à zero. A notação de simetria leva em con¬ 
sideração a simetria do estado molecular. A tabela 2.1 apresenta as três notações 
para alguns orbitais moleculares de moléculas diatômicas, onde cada letra grega 

. |m| 0 12 3 ... 

representa o valor absoluto de m, ou sei a: 

letra a n 5 <fi ... 

As letras s, p, d, f, g, ... representam os valores de /, responsável pela paridade 
do orbital da molécula de conforme a tabela (2.2). 

Os subscritos g e u, designam a paridade da molécula. Se forem pares ou ímpares 
com respeito à simetria de inversão da molécula são denotados pelos símbolos g 
(gerade ) e u (ungerade) 1 , respectivamente. Para que esta classificação seja válida, é 
necessário que o sistema contenha um centro de inversão. Por este motivo a paridade 
somente se aplica a moléculas diatômicas homonuclcares. Então, os estados molecu¬ 
lares relacionados aos estados atômicos s, d, g... possuem simetria g, enquanto que 
os relacionados aos estados p, f, h... possuem simetria u. O simbolo * é usado para 
representar estado anti-ligante ou a ausência dele para representar estado ligante. Se 
existe um estado ligante, existe um valor da distância internuclear R para o qual a 
energia eletrônica é mínima e a configuração molecular permanece estável, também 
chamado de orbital ligante. Quando a energia eletrônica não atinge um mínimo para 


1 Gerade e ungerade são palavras alemãs que significam par e ímpar, respectivamente. 
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quaisquer valores de R, resultando na dissociação da molécula, chama-se de orbital 
anti-ligante. 

2.2.4 Diagrama de Correlação 

O problema consiste em construir um diagrama de correlação que conecta o átomo 
unido (R=0) e átomos separados (R = 00) com os níveis de energia dos orbitais 
moleculares (OM). As propriedades da função de onda dependem da separabilidade 
da equação de Schrõdinger em coordenadas elípticas, pois elas podem ser especifi¬ 
cadas pelo número de nodos em cada coordenada, os números quânticos nodais rii, n 2 
e n 3. Morse & Stueckelberg (1929) descreveram em seu artigo como os três números 
quânticos surgiram no problema de dois centros: n 1, n 2 e n 3 são os números de 
superfícies nodais correspondentes às raízes da solução (equação 2.67) da equação 
de Schrõdinger (equação 2.30) em coordenadas elípticas (A,/q </>), onde ri\ = n\, 
n 2 = riu e n 3 = n 4>- Devido ao fato que o número de superfícies nodais é conser¬ 
vado para todas as distâncias internucleares, esse números quânticos são ”bons”, e 
o diagrama de correlação resultante é denominado diabático, em vez de adiabático, 
que possui regras diferentes. Os números quânticos nodais estão relacionados aos 
números quânticos atômicos da seguinte forma: 


n\ — n\ — n — 1 — 1 

(2.69) 

n 2 = rifj, = í — \m\ 

(2.70) 

n 3 = n 0 = \m\ 

(2.71) 

no limite do átomo unido (Morse & Stueckelberg, 1929). 

No limite de átomos separados os números quânticos 
n \, n' 2 e m', que satisfazem a equação: 

são usualmente dados por 

n s = n[ + n' 2 + \m'\ + 1 

(2.72) 

Como n\ = n\, obtemos: 

n s = n\ + n' 2 + \m'\ + 1 

(2.73) 

onde 



2n ' 2 , gerade 

2 n' 2 + 1, ungerade 
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e 


ricj, = | m! 


(2.75) 


A forma mais geral de construir o diagrama de correlação para o consiste 

em: 


n, = 


n 


n — 


H) 

2 

(l—|m|+l) 


para (£ — \m\) par 
para {£ — \m\ + 1) ímpar 


(2.76) 


A figura (2.3) mostra os diagramas de correlação para os estados n s = 1, n s = 2 


e n s = 3 do Hj. 


átomo unido molécula esferoidal efeito Stark átomos 

separados 


n u n.l.m 

V rym 

n r "1 m 

"s 

6-6 5 0 CJi„ 

0 5 0 CJi<r u 




4 1 

3j/i _ 

03 1 

— 5 tj 

•10 

3 0 

:./« 

0 4 0 

13 0 

,J Ü°o 

5/ tf. 


, -1 3 1 

4/a ~ 

0 1 2 

4/í. 




. 13 1 

-4/,— 

-0 2 1 — 

*/ n 




{- -130 

4/„ 

030 

4 /*. 




\\ 





^200 


-1 2 1 

4r/ • 

— 111 — 

—Ari n a ——_ 








V — 

^10 1 


4 2 0 

-V/„ -— 

— 12 0 — 

—4r/ tf„ 









0 0 2 


4 10 


1 1 0 

tf u 



1 





v ' 








'110 


,•300 

3< — 

2 00 

3i 









v 0 1 1 


3 11 

- 

1 0 1 










' 0 2 (1 


:( 3 2 2 

3da 

00 2 

ü/ í>. 





3 2 I -iW,— —0 I I - M 




110 




3 10 



:%> a ° — 

^ 1 ^ ^ 100 \ 


2 0 0* 


1 0 0 

>111 1 

- 2a tf, — “ 

n n 1 — 

2 t.[i 

2 10- 

- 2/),- 

-u u 1 - 

_(J | Q _ 

-2/) ff D 

010 


2 /*i 


- 2;j tf. ----- 

1=—000 — 

— la 

1 00 

- Is — 

, , , 

- la tf, —- — 




Figura 2.3: Diagrama de correlação em detalhe para os estados n s — 1, n s — 2 e 
n s = 3 do [Adaptado de Allard et al. (1998b)] 
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2.2.5 Análise numérica 

Cada raiz de qualquer uma das relações de recorrência, (2.48) ou (2.61), especifica 
uma relação funcional entre A e p. A solução comum de um par de tais raízes, uma 
de cada relação, resulta nos autovalores de A e p para um certo estado quântico 
(Wallis & Hulburt, 1954). Portanto, o método numérico consiste em resolver o 
sistema de equações não lineares fi(A,e,R) e f 2 (A,e,R) a partir do método de 
Newton-Raphson para sistemas de equações não lineares, implementado em lin¬ 
guagem FORTRAN no nosso caso. 

A estrutura básica e sequência do nosso programa implementado em FORTRAN 
é a seguinte: 

1. Inicialização dos valores de n,l,m,£, A, e R; 

2. Incremento da distância internuclear R; 

3. Gerar as matrizes e efetuar o cálculo dos determinantes obtendo as funções 

MA,£,R)ef 2 (A,£,R); 

4. Resolução do sistema de equações não lineares via método de Newton-Raphson. 

Os valores iniciais dos estados moleculares do são obtidos por aproximação 
em função de R: Quando R —» 0 a molécula comporta-se como um átomo hidrogenoide 
de número atômico Z = 2 com energia eletrônica dada por e = —Aj e A = 1(1 + 1). 

A outra aproximação é quando R —* 00 , e nesse caso a molécula comporta-se como 
um átomo hidrogenoide com número atômico Z — 1 e energia eletrônica dada por 
£ = ~ 2 ^ 2 - Para uma varredura completa de todas as possíveis linhas satélites, pre¬ 
cisão nos cálculos e convergência da solução, o incremento em R adotado em nosso 
cálculo foi de 0,01 u.a., a partir de um valor inicial de R = 0, 01a 0 , onde a 0 é o raio 
de Bohr. As matrizes criadas possuem ordem 25 devido ao critério de convergência 
do programa (IO -21 ) que adotamos. Os valores obtidos para a energia potencial do 
H 2 , possuem discrepâncias somente na última casa decimal (ICC 13 ) quando com¬ 
parados aos valores obtidos por Madsen & Peek (1971); nosso resultado deve ser 
mais preciso. A tabela (2.3) mostra alguns valores do potencial E, (E — e + A), em 
função de R, publicados por Madsen & Peek (1971) comparados com os resultados 
obtidos nesse trabalho para o estado fundamental do H^- 

A figura (2.4) mostra as curvas de energia eletrônica com ( E ) e sem (e) repulsão 
nuclear do estado fundamental (lscr s ) do a partir dos dados obtidos nesse tra¬ 
balho. Em R y 00 , a repulsão nuclear é0eE = e = — \E h . 
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Tabela 2.3: Energia potencial para o estado fundamental do Il.j 


R(au) 


E 


Mareio 


E 


Madsenl971 


I, 0 
2,0 

II, 0 
20,0 
40,0 
90,0 
111,0 


-0,451786313378449655431 

-0,602634214494942477280 

-0,500299229370435249997 

-0,500014259330904820011 

-0,500000881084767589968 

-0,500000034308793139814 

-0,500000014825746297298 


-0,4517863133781 

-0,6026342144949 

-0,5002992293704 

-0,5000142593309 

-0,5000008810848 

-0,5000000343088 



Figura 2.4: Energia eletrônica com ( E ) e sem (T) repulsão nuclear do estado funda¬ 
mental (lscx 5 ) do a partir dos dados obtidos. 

A curva E tem um mínimo em R = 2, OOao = 1,06 Ã, indicando que o estado 
eletrônico fundamental do H^“ é um estado ligante estável. O valor obtido para e 
a 2,OOa 0 foi de -1,102634214494942144213 E h e somando o termo devido a repulsão 
nuclear -b, obtém-se E = —0,602634214494942477280 E h , comparável com —\E h 
para R —* oo. Consequentemente, a energia de dissociação do estado fundamental 
é D=0,10263421 E h = 2,79282056 eC = 64,40395570 kcal/mol, ou seja, 17,02% da 
energia total ( E ) na separação de equilíbrio. 

As figuras (2.5) e (2.6) mostram a energia eletrônica e a energia potencial re- 
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spectivamente dos estados com n s = 1 do H^. 



As figuras (2.7) e (2.8) mostram a energia eletrônica e a energia potencial re¬ 
spectivamente dos estados com n s = 2 do H^". 





Energia (E h ) Energia (E h ) 
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Figura 2.7: Energia eletrônica (e) para os estados com n s = 2 do Hij - . 
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Figura 2.8: Energia potencial (E) para os estados com n s = 2 do . 








Capítulo 3 
Dipolo 


Os momentos de dipolo elétrico produzem alterações na intensidade dos satélites nas 
asas distantes dos perfis de linha (Allard et ah, 1999). Os perfis de linha mostram 
que a amplitude dos satélites depende da intensidade do momento de dipolo elétrico 
na distância internuclear onde o satélite é formado. Como consequência, a variação 
do momento de dipolo elétrico deve ser levada em conta para se obter resultados 
confiáveis, pois são usados como diagnósticos de parâmetros estelares e de plasmas, 
portanto, não podem ser negligenciados. 

3.1 Momento de dipolo elétrico do íon m 

O momento de dipolo elétrico para transições entre estados eletrônicos ligados do íon 
molecular Ht para uma dada separação internuclear pode ser obtido analiticamente, 
usando as soluções de Hylleraas (1931) da função de onda. 

3.2 Funções de onda do m 

Escrevendo a equação de Schrõdinger para o íon molecular Ht em termos de coorde¬ 
nadas elípticas confocais (A,/q 0), obtém-se a equação (2.30), que é completamente 
separável e cuja solução pode ser escrita como um produto de três funções: 

tt(A,/z,0) = A(A)M(/z)$(0) (3.1) 

Substituindo a solução (3.1) na equação de Schrõdinger escrita em coordenadas 
elípticas (2.30), obtivemos as equações (2.32, 2.33 e 2.34). 
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A solução da equação 2.32, parte azimutal, é da forma: 

<f>(0) = e (±ím<w , (3.2) 

A solução para a equação 2.33, parte angular, é dada por Hylleraas (1931), 

/ 

aí(m) = (3.3) 

S 

como uma expansão em polinómios associados de Legendre P™ +S (fj,), onde o so¬ 
matório somente se estende sobre os valores de s pares ou ímpares, dependendo 
da paridade do estado e onde os coeficientes c s são obtidos via frações contínuas, 
conforme equação (2.52). 

A parte radial será resolvida conforme Hylleraas (1931) e não mais de acordo 
com Jaffé (1934), como foi feito para obtenção da energia eletrônica e. Portanto a 
solução da equação (2.34) é dada por: 

A(A) = (A 2 - 1 r /2 e- s/2 p(í) (3,4) 

í = 2p(A - 1) (3.5) 

oo 

»<0 = E(STÕi*S«(í) (3 ' 6) 

t=0 ' ’ 

onde LZ+t (£) são os polinómios de Laguerre associados e, devido às suas pro¬ 
priedades, são mais adequados nos cálculos dos momentos de dipolo elétrico. Os 
coeficientes g t são obtidos a partir de uma nova relação de recorrência gerada pelo 


uso da solução de Hylleraas (1931) na parte radial, 

<Xt9t+ 1 - Pt9t + lt9t- 1 =0 (3.7) 

a t = (t — a — m)(t + m + 1) (3.8) 

Pt = {2í 2 + 2t(2p - a) - [A - p 2 + 2pa + (m + 1 ){m + cr)]} (3.9) 

7t = t(t — 1 — cr) (3.10) 
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3.3 Componentes do momento de dipolo elétrico 

Escolhendo, conforme a figura (2.1), a origem do sistema de coordenadas no centro 
do segmento de reta que une os dois prótons e tomando essa reta como eixo z, os 
elementos da matriz a ser calculada dos componentes do momento de dipolo elétrico 
de uma transição entre os estados i e j, D,j , podem ser escritos como: 

Dij = e {i\r\j) = eQij( r) , (3.11) 

onde e representa a carga do elétron, r é o vetor posição do elétron parametrizado em 
função de R (distância internuclcar) e Qíj( r) os elementos da matriz. Os elementos 
da matriz Qíj( r) dos vários componentes do momento de dipolo elétrico, de acordo 
com Ramaker & Peek (1973), são especificados por: 

Qij( r ) = f ^i(A,^»r^(A,/r»dr (3.12) 

onde Ü/j e são as funções eletrônicas na aproximação de Born-Oppenheimer dos 
estados i e j, e dr é o elemento de volume. Qíj( r) é obtido inserindo as equações 
(3.2), (3.3) e (3.4) na equação (3.12). Portanto os elementos da matriz Qij(r) podem 
ser escritos, de acordo com Herman & Wallis (1956), como função das coordenadas 
x, y and z como: 


r»27r 


Qi,j{x) = NíNj 


'-Ui 


R 3 

4/*(À, /i, (j )) {x) 4/j(A, /i, <f>)— (A 2 — y 2 ) d\dfid(j) (3.13) 


Qi,j (*^) 


'•2n 


Ni Ni 


J i 


'k*(A,/i,(/))'k i (A,/i,0) - cos0(A 2 


l)i(l 



x 



(A 2 — y 2 ) dXdydcj) 

(3.14) 


-2tt 


Qi,j(x) = NíNj 


x 


’-iJ i 

,2-í ’ 


-R 4 

^*(A,/^»^(A,/r» — cos0 x 


(A 2 — 1)2(1 — ^n 2 ) 2 A 2 — (A 2 — 1)2(1 — ^n 2 ) 2 


dXdydcp 


(3.15) 
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QiA*) = ^4 





X 


a*(a)m;(/,)$*( 0 )a,.(a)m, (/,)$,( 0) cos</> x 
(A 2 - 1)^(1 - yu 2 )^A 2 - (A 2 - 1)^(1 - /x 2 )^// 2 dXd/idcj) 


0 J-lJl 
2 


(3.16) 


Agrupando os termos, obtemos para Q h j(x): 


Qi,j (*^) 


$*(</>) cos(<£)$, (<£)# x 



1 1 

—, 0, —, 2 1 - O* ? I - 2, 0 

o ’ ’ o ’ / * J \ o ’ ’ o ’ 


(3.17) 


analogamente para Qij(y), 


Qi,j(y) 


&*( 4 >)sen((p)&j( 4 >)d(t) x 



—, 0 , 2 - íí i7 - -, 2,-,0 

O * * O * L iJ \ o * * o * 


(3.18) 


e para Q tú (z), 


Qij(z) = N,N, (jj [’ VA<t>)*,m<t> X 
x[í2 w (0,l,0,3)-n jJ .(0,3,0,l)] 


(3.19) 


onde íljj(a, b, c, d) é definido por: 


Qij(a,b,c,d) — I —/j, 2 ) a iu, b Mj(ii)dn x 


1 


x / A*(A)(A — 1) C A Aj(A)dA 


(3.20) 


As constantes de normalização podem ser calculadas como: 


Nr = tt [%( 0, 0,0, 2) ~íi M (0, 2, 0, 0)] x /$*(</>) {(j>)d(j> (3.21) 
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fe(0,0,0,2) 


ÍX,-j(0, 2, 0, 0)] x 


<f>*(0)%(0)d0 


(3.22) 


Para a solução completa dos valores de Qíj( r) faz-se necessário calcular as in¬ 
tegrais que compõem as expressões para Qíj(x), Qij(y ) e Qíj(z ) de cada transição 
Podemos classificar as integrais em três tipos, a saber, as que dependem de 4>, 
as que dependem de fi e as que possuem dependência em À. No que segue trataremos 
dos três tipos e suas respectivas soluções. 

As integrais, 



$*(</>) COS 





2tt 




J 0 

que dependem de 0 e cujos valores são diferentes de zero ao substituir as funções 
$(0) por e im< ^ são: 


f 2n ( x\ ( <x\a I n ’ para m ' = m = 2 ’ 3 ’'' ' 

/ sen{mç))sen{m (p)d(p — < 

/o |0, para os demais casos 


(3.23) 



7T, para m' — m — 1, 2, 3, • • • 1 

/ cos(ra</>) cos (m'(j))d(f) = < 

27t, para nn! — m — 0 / 

/o 

0 , para os demais casos J 


(3.24) 


r» 27 r 


sen{m4>)sen{m!'</>) cos <j)d(j) = 


|, para m! = m + l,m^0 
0, para os demais casos 


(3.25) 




Dipolo 


46 


r*27r 


cos (m<t>) cos (m'4>) cos 0d0 = 


para m! — m + 1, m ^ 0 
7T, para m' = m + 1, m = 0 
0, para os demais casos 


(3.26) 


-2tt 


sen{m<j )) cos {m!4>)sen4>d4> = 


para m ; = m + l,m/0 
0, para os demais casos 


(3.27) 


n2n 

para m' — m + 1 , m ^ 0 1 

/ cos(m0)5en(m'çí))5en0dçí) = < 

7T, para m' = m + 1 , m = 0 / 

/ 0 

0 , para os demais casos J 


(3.28) 


As integrais dependentes de (/x) são: 


JÇW(> 


'-1 


H 2 ) a H b Mj{n)dn = A 


cuja solução depende dos estados (i,j) e dos valores de (a, b), da seguinte forma: 
Para i=jea = b = 0: 




2(2 m + s)! 

(2 m + 2s + l)s! 


Para i = j, a = 0 e 6 = 2: 




2(2m + s)! 
(2m + 2s + l)s! 


0,1 (c«—2 d" C s 


i a 2,m,s "P C s _)_2,ifl2, 


m,s+2/ 


Observando as soluções das integrais dependentes de 0, equações 3.23 a 3.28, 
somente os valores em que m' = m e m' = m + 1 são não nulos. Os estados i e j 
são escolhidos de tal maneira que m! > m. 

Para o caso em que m' — m + 1, 
com i 7 ^ j, a — \ e b — 0: 




2(2 m + s)! 
(2m + 2s + l)s! 


^m,s—2^s—2,j ) 
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e com i ^ j, a = \ e b — 2: 
2(2 m + s)! 


* = E 


|_(2m + 2s + l)s! 


ds,i 0 , 2,m,s—2^rn, S — 4*-"S—4 jd" (®2,m,s—2^m,s—2 ^‘2,m,s^m,s—2)^s—2,_ 


^~( a 2 ,m,s^m,s a 2,m,s+2^m,s) C s,j d~ ( a 2,m,s+2^m,s+2 c s+2 j] 

Para o caso em que m! = m, com i ^ j, a — 0 e b — 1: 

2(2 m + s)! 


* = E 


(2m + 2s + l)s! 
e com * ^ j, a — 0 e b — 3: 

2(2m + s)! 


Cs,i (®l,m,s—l^s—lj “h ffi s-f-iCs+ij) 


A = E 


(2m + 2s + l)s! 


c s ,i(a3 ,m,s— 3 c s-3,j+ a 3 ,m,s— l^s -lJ+°3 ,m,s +l c s+l,j _ l _a 3,m,s+3 C s+3j) 


Os termos ai_». 3 )m)S e b myS são os coeficientes das relações de recorrência dos 
polinómios associados de Legendre (Abramowitz & Stegun (1972), Arfken & Weber 
(1995)): 


A^tti+sÍaO ®l|»*iS^ra+s+l(^) ^ 


^m+s (aO ~ ®2,m,s-fm+s+2 (aO A 


a' p m (,,) n " p< 

^2 777, ,S 7 T},A-fí\r J ) ^9 m «-* m-Uc- 






(i - ^PZtLM = V.,-P” + ,(í») - fC.CW», 


e são dados por: 


(s + 1) 

1)m ’ s (2m + 2s + 1) ’ 

, _ (2 m + s) 

0l ’ m ’ s “ (2m + 2s + 1) ’ 
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_ (s + 2) 

2 ’ m ’ s ~ (2m + 2s + 3)' 1>m ’ s ’ 

/ (2m + s + l) , (s) , 

a2 ’ m ’ s ~ (2 m + 2s + 3) 1>m ’ a + (2m + 2s - 

„ _ (2 m + s - 1) , 

° 2 ’ m ’ s " (2m + 2s — 1) 0l ’ m ’ s ’ 

_ ( s + 3) 

3 ’ m ’ s_ (2m + 2s + 5) 2 ’ m ’ s ’ 

z (2 m + s + 2) | (s +1) , 

a3 ’ m ’ s = (2m + 2s + 5) a2 ' m ' s + (2m + 2s + l) a2 ’ m ’ s ’ 

// = (2m + s) , (f ~ 1) „ 

3 ’ m ’ s (2 m + 2s + 1) 2 ’ m ’ s (2 m + 2s - 3) 2 ’ m ’ s ’ 

_/// = |2m + s-2[ „ 

3 ’ m ’ s (2m + 2s - 3) 2 ' m ' s ’ 

(2 m + s + 1) (2 m + s + 2) 

6m ’ s = (2m + 2s + 3) ’ 

, = (s + l)(s + 2) 

m ’ s (2m + 2s + 3)' 

As integrais dependentes de À são: 


A*(À)(A 2 — l) c X d Aj(X)d\ = E 


onde a solução depende dos estados (i,j) e dos valores de (a, 6), da seguinte forma: 
Para i = j e a = b = 0: 


^ ^ ^ 9t,i9t’,i ^ ' 

t,f v=o 

Para i = j, a = 0 e b = 2: 



0 + m,pi,pj) 




í,t' 


i >=0 


X 


+ m + 2,pi,pj) + 2/™; m (u + m + 1,Pí,Pj) + + m,p i ,p j )_ 
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Para i ^ j, a — \ e b = 2: 

oo m +1 / . \ 

) 2 ” _ ” +lx 

t,t' ?>=o v 7 

x [r t y n+1 (v + m + 3,pi,pj) + 2/™; m+1 (u + m + 2,pi,pj) + /™; m+1 (v + m + l,Pi,Pj)\ 
Para i 7 ^ j, a — \ e b — 0: 

00 m +1 , 1 1 \ 

s 9t,i9t>,% ( m v ) 2 m ~ v+1 I™r +1 (v + m + 1, pi, Pj) 

t,t' v=0 x V 7 

Para i 7 ^ j, a — 0 e b — 3: 


x 


==e^,.e 

í,t' v=0 V 7 

x [J™;' n (u+m+3, Pi, pj)+3r t y n (v+m+2, p h pj)+3I™ i ’, m (v+m+l,p i , p j )+I^ rn (v+m, p h pj)\ 
Para i ^ j, a — 0 e b — 1: 


-E—EQ 2 

t,t' u =0 ' 7 

x [/™; m (v + m + 1, Pi, Pj) + i$> + m, Pi, Pj)} 

Nas expressões acima os termos (“) são números binomiais e os valores de pt e 
Pj são obtidos a partir de p 2 = — A integral J )";' 77 ( q,Pi,Pj ) pode ser escrita 
função de p = À — 1 = A da seguinte forma: 


x 


em 


jm.m / \ 

h,t' (e> Pú Pj) = 


roo í -| \m-\-m 

/ 1 > _ n Q p -ÍPi+Pj)P x 

/ 0 (m + t)\(m'+ t')\' 

x L™ +t (2 pip)L™, +t , (2 pjp). 


(3.29) 


A solução da integral (3.29), de acordo com Herman & Rubin (1955), possui a 
forma: 


t v 


Çt' n \(LPi,Pj) = ( Pi+Pj ) (1+g) EE(-i) 




0 i ;=0 


X 


t + m\ /P + rn! 


t — u 7 V P — v 


x 


T(1 + ç + tt + v)(2pi) u (2pj) 1 
u\v\(pi + Pj) íl " H ' 


(3.30) 
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Tendo calculado os valores de Qíj( r), obtemos o valor de D tJ da seguinte maneira: 


Dij | 2 — e [\Qij(x )\ 2 + | Qi,j(y )\ 2 + | Qi,j(z )\ 2 


(3.31) 


3.4 Análise numérica 

As regras de seleção dipolar são: 

1. Transição entre estados com paridade oposta; 

2. A variação em \m\ deve ser zero ou unitária (|m — m'\ = 0,1); 

3. Não existe restrições na mudança em n u . 

Baseado nas regras de seleção, foram calculados os valores do momento de dipolo 
elétrico, equação (3.31), para as transições das séries de Lyman e Balmer do até 
n = 10, com separações nucleares até 300 a 0 . As convergências das séries envolvidas 
no programa escrito em linguagem FORTRAN foram distintas, ou seja, cada série 
teve que ser analisada individualmente até a sua convergência. Para as equações 
3.3 e 3.4 serem soluções apropriadas das funções de onda as séries devem ir ao 
infinito. Dentro da precisão numérica proposta, a convergência ocorreu para s = 31 
e t — 10. A tabela 3.1 mostra os valores dos momentos de dipolo elétrico desse 
trabalho comparados com os obtidos por Ramaker & Peek (1973) para algumas 
transições, demonstrando que os cálculos descritos aqui apresentam grande precisão 
e cobrem uma vasta distância internuclear. 

A figura 3.1 mostra as transições listadas na tabela 3.1 dos momentos de dipolo 
elétrico calculados nesse trabalho. 
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Tabela 3.1: Tabela de comparação entre os dados de Ramaker & Peek (1973)! e 
esse trabalho^ dos momentos de dipolo elétrico para Lyman-a do H^, em unidades 
atômicas. _ 


R{a o) 

Transição 

Dij(R)' 

DiARf 

1,0 

2pn u -lsa g 

0,5494 

0,5493774 

11,0 

2jm u -lsa g 

0,7169 

0,7168555 

60,0 

2jm u -lsa g 

0,7450 

0,7449931 

1,0 

3pa u -lsa g 

0,1651 

0,1651107 

11,0 

3pa u -lsa g 

0,4859 

0,4858565 

80,0 

3pa u -lsa g 

0,5207 

0,5207119 

1,0 

4fa u -lsa g 

0,001055 

0,0010541 

11,0 

Afa u -ls(jg 

0,6609 

0,6608712 

80,0 

Afa u -lsa g 

0,5326 

0,5325508 

1,0 

2sa g -2pa u 

1,241 

1,2406089 

11,0 

2sa g -2pa u 

0,4844 

0,4844248 

80,0 

2so g -2pa u 

0,5207 

0,5207119 

1,0 

3da g -2pa u 

1,010 

1,0104411 

11,0 

3da g -2pa u 

0,6349 

0,6349347 

80,0 

3da g -2pa u 

0,5326 

0,5325508 

1,0 

3dn g -2po u 

0,8584 

0,8584398 

11,0 

3dn g -2po u 

0,7935 

0,7935359 

60,0 

3dn g -2pa u 

0,7450 

0,7449931 




D(R)(ea ) 
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Figura 3.1: Momento de dipolo elétrico das transições de Lyman-a:, listadas na 
tabela 3.1, devido ao átomo de H perturbado por colisões com prótons. O D(R ) é 
dado aqui em unidades atômicas e a distância internuclear R em a 0 (raio de Bohr). 










Capítulo 4 

Perfis e satélites 


4.1 Análise teórica 

Na introdução já foi discutido que a técnica espectroscópica fornece medidas precisas 
dos parâmetros atmosféricos estelares, T e s e log g, permitindo dessa forma a com¬ 
paração relativa desses parâmetros entre estrelas distintas. No entanto, os valores 
absolutos dos parâmetros atmosféricos podem ser imprecisos devido às incertezas na 
física incluída nos cálculos do modelo (e.g., Tremblay & Bergeron, 2009). A principal 
fonte de alargamento, na maioria das estrelas anãs brancas DA, é devido à interação 
das partículas com o campo elétrico, ou seja, o alargamento Stark. A teoria que 
descreve esse perfil de linha até hoje é a teoria unificada de alargamento Stark de 
Viciai, Cooper e Smith (Vidal et ah, 1970) conhecida como VCS. Essa teoria foi us¬ 
ada desde o início na análise de grandes amostras de estrelas anãs brancas DA e é a 
base da maioria dos resultados aceitos das propriedades globais dessas estrelas (e.g., 
Tremblay & Bergeron, 2009). Lemke (1997) disponibilizou novas tabelas do alarga¬ 
mento Stark para as séries de Lyman, Balmer, Paschen e Brackett do hidrogênio, 
calculadas a partir da teoria de VCS. Essas novas tabelas cobriram todas as possíveis 
transições, temperaturas e densidades eletrônicas encontradas nas atmosferas de es¬ 
trelas anãs brancas DA. Antes do trabalho de Lemke (1997) os perfis de linha eram 
obtidos por extrapolação para altas densidades ou utilizando teorias de alargamento 
menos precisas, como de Edmonds et al. (1967), para linhas superiores a He da 
série de Balmer. Mas para temperaturas da ordem de 10 000K para as linhas de 
Balmer e abaixo de 30 000K para as linhas de Lyman, outros tipos de alargamento 
de linha devem ser incluídos no modelo de atmosferas. Uma das aproximações mais 
expressivas tem sido o trabalho de Allard et al. (2004) sobre opacidade de linhas 
quase moleculares. 
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Existem diferentes aproximações para o cálculo dos perfis de linha; as principais 
são: 


• aproximação quase-estática; 

• aproximação quântica, 

• teoria unificada. 

Na aproximação quase-estática, o coeficiente de absorção, K[y) (medido em 
cm -1 ), causado pelas absorções induzidas por colisões (CIA) 1 (Margenau & Lewis, 
1959, Allard & Kiclkopf, 1982, Rohrmann et ah, 2011), pode ser escrito como: 


K(v) = n Hi crtot(v) 

Onde cr íoí (i/) é o coeficiente de absorção total, dado por: 


(4.1) 


tftoíH = y (4.2) 

ij 

Os valores são os pesos de cada transição de um estado i para um estado j de um 
átomo ou molécula e (7 l3 (v) é a seção de choque da radiação absorvida na respectiva 
transição. Para o caso de H — H + , o coeficiente é dado por (Bates, 1951, 1952, 
Rohrmann et ah, 2011): 


a 


o 


W) = 


16 , 

—n n p a 


hu 


\dVij/dR | 


í) IJ R~e 


.2 Uixn p R 3 /S) -pEi(R) 


(4.3) 


onde f3 = -jjf, k é a constante de Boltzmann, e T é temperatura do gás. O n p é 
a densidade média de perturbadores no gás, V l3 é a diferença da energia potencial 
associada com uma transição com uma energia hv = V tJ = V 3 (R) — Vi(R) para uma 
separação internuclear R, D í3 (R ) é o momento de dipolo, e Ej(R ) corresponde a 
mudança na energia (com relação a R = oo) no estado i do H^. Quando existe 
um extremo na curva correspondendo à diferença nas energias potenciais entre um 
estado superior e inferior,^- é nula. Portanto, se a derivada é nula, a expressão do 
perfil é divergente, levando aos satélites nas asas dos perfis. 

Na aproximação quântica, a seção de choque para processos associados espontâneos 
e estimulados é dada por (Zygelman & Dalgarno, 1990), 


1 Collisional Induced Absorptions 
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< E ) = Y.Y, 


64 7T 5 z/ 3 


X 


3 c 3 K 2 1 — exp(—hv/kT) 
p[N + (N + l)M\ N , +1 . k N ,) 


AT / H 

N v 
r2 


(4.4) 


onde E é a energia relativa da colisão, k é o número de onda do movimento relativo, 
pê a probabilidade de aproximação no estado eletrônico inicial, N' ê o número 
quântico rotacional inicial, v" é número quântico vibracional final, e M é o elemento 
da matriz de dipolo elétrico. 

A teoria unificada foi desenvolvida em Allard et al. (1999), e uma discussão detal¬ 
hada é apresentada nesse artigo. A expressão fundamental do espectro normalizado 
F v {Au) é calculada a partir da transformada de Fourier (Allard et ah, 2009). 


-| r+oo 

F v (Au) = -Re / $(s)e- lAus ds. (4.5) 

71 J 0 

onde <F(s) é a função de autocorrelação dipolar (Allard et al., 1999). 



Figura 4.1: Comparação entre as três aproximações do perfil de linha de Lyman-a: 
normalizado: aproximação quase-estática (linha tracejada), aproximação quântica 
(linha pontilhada) e a teoria unificada (linha sólida) [Adaptado de Allard et al. 
(1999)] 


A figura 4.1 mostra a comparação entre as três aproximações, mencionadas 
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acima, do perfil de linha de Lya como calculado por Allard et al. (1999). A 
linha sólida representa a aproximação da teoria unificada, que possui uma resolução 
determinada pela transformada de Fourier (equação 4.5), enquanto que a aprox¬ 
imação quântica (linha pontilhada) e a aproximação quase-estática (linha tracejada) 
fornecem valores semelhantes, mas perfis distintos devido ao cálculo da derivada pela 
diferença numérica na aproximação quase-estática, e por expansão em série na aprox¬ 
imação quântica. Os três perfis são semelhantes, mas a aproximação quase-estática 
melhor representa a forma correta global das asas da linha e a posição dos satélites, 
em comparação com dados publicados anteriormente (Allard et ah, 1999). 

Como os objetivos principais desse trabalho são os cálculos dos momentos de 
clipolo elétrico, os potenciais do íon H^", os cálculos da forma global dos perfis de 
linha e a localização dos satélites, usamos a aproximação quase-estática para calcular 
os perfis de linha. Infclizmente a aproximação quase-estática não fornece uma boa 
representação do centro das linhas dominado pelo efeito Stark. Portanto, há uma 
necessidade de uma melhor modelagem das partes centrais de cada transição. 

Nosso interesse é encontrar a posição correta dos satélites nas asas do perfil de 
linha e sua forma geral, ou seja, as divergências no perfil, e levar esses resultados 
para a atenção dos astrofísicos. Se eles são observados, são devidos a prótons. 

4.2 Aproximação quase-estática 

A abordagem teórica usada nesse trabalho é baseada na absorção quase-molecular H 
- H + calculada na aproximação quase-estática, da mesma forma que Rohrmann et al. 
(2011). O método de cálculo segue o procedimento de Rohrmann et al. (2011), para 
perfis de linha das séries de Lyman e Balmer em uma aproximação quase-estática: 

• Primeiramente identifica-se cada par de níveis (i,j) da quase-molécula H - H + 
que contribui para a formação de uma dada linha. Somente transições entre 
estados u ( ungerade ) 2 e g ( gerade ) 3 são permitidas (transições dipolares), isto 
é, levam a emissão de um fóton, ou seja, com paridade oposta e possuindo 
| m — m! |= 0 ou 1 como clipolo permitido. 

• O segundo passo é o calculo da seção de choque efetiva equação (4.3) 

usando os potenciais e momentos de clipolo elétrico para as transições já iden¬ 
tificadas. 

2 Palavra alemã para ímpar 
3 Palavra alemã para par 
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• Para o cálculo correto de ay,-(V) precisa-se analisar as frequências v que surgem 
a partir da diferença da energia hv = Vj — Vi entre os níveis que contribuem 
para a formação da linha. Deve-se levar em conta que distâncias internucleares 
distintas podem produzir frequências iguais ou extremamente próximas. 

• A seção de choque efetiva é a soma de cada uma das contribuições de todas as 
possíveis transições. 

• O perfil de linha é a soma das seções de choque efetivas <Jij(v) multiplicadas 
pelo peso estatístico de cada uma das respectivas transições permitidas. 

O perfil calculado é então convoluído com uma distribuição Maxwelliana de ve¬ 
locidades para levar em conta o efeito do alargamento térmico Doppler. Os perfis 
de linha e satélites das séries de Lyman e Balmer foram calculados para densidades 
únicas de 10 16 até 10 18 cm -3 típicas das atmosferas de estrelas anãs brancas e que 
podem ser produzidas em laboratório em experimentos Z-pinch. A última etapa 
para a obtenção do perfil real é somar o perfil obtido ao perfil devido somente ao 
alargamento por elétrons, alargamento Stark, o qual não possui satélites. 

4.3 Série de Lyman e Balmer 

As linhas de Balmer do hidrogênio atômico são proeminentes na maioria dos espec¬ 
tros estelares óticos bem como as linhas de Lyman na maioria das estrelas quentes o 
suficiente para terem um fluxo de energia no ultravioleta observável. A interpretação 
correta dos espectros e, por consequência, a determinação dos parâmetros estelares, 
necessita do conhecimento da forma do coeficiente de absorção da linha. O objetivo 
dos cálculos apresentados aqui é fornecer perfis de linha para as condições físicas 
nas atmosferas de estrelas anãs brancas, plasmas e experimentos do tipo Z-pinch. 
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4.3.1 Lyman-o 

O perfil total de Lyman-o; depende de seis transições individuais da molécula de H^, 
todas listadas na tabela 4.1. A tabela 4.1 também exibe o comprimento de onda 
das linhas satélites, bem como o nível superior e inferior da transição que gera o 
respectivo satélite e a distância entre H — H + para a qual o extremo na diferença 
de potencial ocorre. 


Tabela 4.1: Transições permitidas e satélites devidos a colisões de H-H + para Lyman- 
a na distância R(Â). 


N° 

Nível superior 

Nível inferior 

Ais h + (cm x ) 

^h+ (Â) 

RlA) 

1 

2pTT u 

lsa g 

-1628,54 

1 240,56 

4,72 

2 

3pa u 

lso g 

- 

- 

- 

3 

4 fo u 

lSOg 

-1172,40 

1233.58 

11,07 

4 

2sa g 

2 pcr u 

8139,91 

1106,47 

2,65 

5 

3da g 

2 pc>u 

-11055,90 

1404.87 

4,56 




16 001,33 

1017.93 

1,39 

6 

3d7T g 

2 pcr u 

21384,28 

965,05 

1,67 


A figura 4.2 exibe as seis transições listadas na tabela 4.1. Nessa figura é possível 
observar as distâncias onde ocorrem os extremos da diferença de potencial, que são 
responsáveis pela criação dos satélites no perfil de linha. 

Os momentos de dipolo elétrico de cada uma das transições que contribuem para 
Lyman a, são exibidos na figura 4.3. 

De acordo com o método de Rohrmann et al. (2011), para obtenção do perfil de 
linha primeiramente calcula-se a contribuição de cada transição, as quais podem ser 
observadas na figura 4.4 para uma temperatura de 10 000 K, com um número de 
perturbadores n p = 1,0 x 10 17 cm -3 de H + . O indicador de que cada contribuição 
está correta é a linha de corte que cada uma delas apresenta. A figura 4.5 mostra 
todas as transições exibidas na figura 4.4, mostrando como cada uma delas contribui 
na formação do perfil de linha. 
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Figura 4.2: Diferença de potencial das transições permitidas em função da distância 
interatômica da molécula de H^" de Lyman-a. 
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Figura 4.3: Momento de dipolo elétrico das transições de Lyman-a devido ao átomo 
de H perturbado por colisões com prótons. O D(R ) é dado aqui em unidades 
atômicas. 
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Figura 4.4: Transições que contribuem para a formaçao de Lyman-a mostradas 
individualmente para T=10 000 Ken p = l,0x 10 17 cm“ 3 
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Figura 4.5: Todas as transições que formam o perfil linha de Lyman-a juntas, ex¬ 
ibindo a contribuição de cada uma delas para T=10 000 K e n p — 1, 0 x 10 17 cm -3 
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l(À) 


Figura 4.6: Comparação entre os perfis de Lyman-a, devido ao alargamento Stark 
e colisional do H — H + para T=10 000 K e n p = 1, 0 x 10 1 ' cm -3 . Note que as asas 
das linhas não são idênticas assintoticamente. 
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X(À) 


Figura 4.7: Comparação entre os perfis de Lyman-a, devido ao alargamento Stark 
e colisional do H — H + convoluído com o efeito Doppler para T=10 000 K e n p = 
1, 0 x 10 17 cm~ 3 
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4.3.2 Lyman-/? 

O perfil de linha de Lyman-/3 depende de dez transições individuais da molécula de 
H^, que estão listadas na tabela 4.2. A tabela 4.2 lista todos os comprimentos de 
onda dessas transições e as possíveis linhas satélites. As linhas satélites de Lyman-/3 
foram observadas nos espectros da estrela anã branca DA Wolf 1346 com o telescópio 
HUT (Hopkins Ultraviolet Telescope ) (Koester et ah, 1996). Koester et al. (1998) 
observaram satélites característicos nos comprimentos de onda 1060Â e 1 080Ã nos 
espectros de quatro outros alvos obtidos com o ORFEUS Orbiting and Retrievable 
Far and Extreme Ultraviolet Spectrometer. Essas linhas satélites do perfil de linha 
de Lyman-/3 são devidas às contribuições das transições (4/7t u -1sct s ) e (5ga g -2pa u ), 
conforme tabela 4.2. Tais linhas satélites também foram calculadas por Allard et al. 
(1998b). A figura 4.8 exibe a diferença entre os potenciais das transições que con¬ 
tribuem para a formação do perfil de linha de Lyma n-/3, listadas na tabela 4.2. 
Na figura 4.9 são graficados os momentos de dipolo elétrico de cada uma das dez 
transições. Observa-se nas figuras 4.8 e 4.9 que quando a diferença de potencial 
produz um extremo e o momento de dipolo elétrico possui declividade positiva, há 
formação de uma linha satélite no perfil de linha. O perfil de linha de Lyman -(3 
calculado para uma temperatura de 10 000 K e número de perturbadores (prótons) 
n p = 1, 0 x 10 17 cm -3 é exibido na figura 4.10. 

Tabela 4.2: Transições permitidas e satélites de Lyman-/3 devido a colisões do átomo 
de H com H + na distância internuclear R(Ã) na qual um extremo no potencial ocorre. 


N° 

Nível superior 

Nível inferior 

Au h + (cm x ) 

^ h+ (Â) 

R( A) 

1 

3p7T u 

lsa g 

- 

- 

- 

2 

4 pcr u 

lSdg 

- 

- 

- 

3 

4fTT u 

\S(Jg 

-3 361,25 

1062,64 

9,85 

4 

5 fa u 

lSGg 

- 

- 

- 

5 

6ha u 

ISOg 

-1037,33 

1037,03 

21,45 

6 

3 SGg 

2 pcr u 

22 200,04 

835,65 

1,96 

7 

Ad(Jg 

2 pc>u 

-633,34 

1032,71 

9,46 




25 868,59 

810,80 

1,55 

8 

Adllg 

2 po u 

28 666,11 

792,82 

1,70 

9 

3go g 

2 pcr u 

-4 899,9272 

1080,31 

12,65 




38 880,59 

733,42 

1,69 

10 

3gn g 

2 po u 

-516,61 

1031,46 

18,89 




38 914,86 

733,24 

1,69 
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Figura 4.8: Diferença de potencial das transições permitidas em função da distância 
interatômica da molécula de de Lyman-/3. 
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Figura 4.9: Momentos de dipolo elétrico das transições que contribuem para 
formação de Lyman-d devido ao átomo de H perturbado por colisões com prótons. 
O momento de dipolo elétrico D(R) é dado aqui em unidades atômicas. 
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Figura 4.10: Perfil de linha de Lyman ~/3 calculado para a temperatura T=10 000 K 
e para uma densidade de perturbadores n p = 1, 0 x 10 17 cm -3 . 
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4.3.3 Lyman -7 

O perfil de linha para Lyman -7 é formado por 14 transições permitidas do íon molec¬ 
ular H 2 devido às colisões de H-H + . Essas transições são listadas na Tabela 4.3, 
bem como a posição e a distância internuclear entre H e H + dos possíveis satélites 
gerados. Um dos satélites predito para ocorrer, em À = 995 Ã, foi observado em 
um espectro do HUT (Hopkins Ultraviolet Telescope) e mais tarde em um espectro 
obtido com ORFEUS (Orbiting and Retrievable Far and Extreme Ultraviolet Spec- 
trometer). Esse satélite em Lyman -7 também foi observado nos espectros das anãs 
brancas CD — 38°10980 (Wolff et al., 2001) e Sírius B (Holberg et al., 2003). Allard 
et al. (2004) calcularam o perfil de linha de Lyman -7 utilizando os potenciais do H 2 
calculados por Madsen & Peek (1971) e assumiram que o momento de dipolo era 
constante para cada transição. Allard et al. (2009) recalcularam o perfil de linha 
para Lyman -7 com a inclusão dos momentos de dipolo elétrico para cada transição, 
mas os dados de dipolo não foram publicados. 


Tabela 4.3: Transições permitidas e satélites devido às colisões H-H + de Lyman -7 na 
distância internuclear R na qual um extremo para a diferença de potencial ocorre. 


N° 

Nível superior 

Nível inferior 

Au h + (cm x ) 

^H+ (^) 

;?(Â) 

1 

4pvr n 

lsa g 

- 

- 

- 

2 

5 po u 

lS(Jg 

- 

- 

- 

3 

5 fTTu 

lS(Jg 

-192,78 

974,62 

16,65 

4 

Qfcr u 

lSdg 

- 

- 

- 

5 

6Im u 

ISOg 

-2188,78 

993,96 

20,69 

6 

7 h<y u 

ISOg 

-230,90 

974,98 

29,66 

7 

OO 

£ 

ISOg 

-784,17 

980,27 

36,06 

8 

4 SVg 

2 \pa u 

30 153,49 

752,16 

1,82 

9 

5 d(7g 

2 pcr u 

32 004,30 

741,83 

1,62 

10 

ÒdlTg 

2 pcr u 

33 520,98 

733,58 

1,70 

11 

6 go g 

2 pcr u 

-1395,38 

986,18 

17,79 




39 005,56 

705,21 

1,69 

12 

6.977 

2 po u 

39 025,78 

705,11 

1,69 

13 

hOg 

2 pcr u 

-2 627,07 

998,31 

25,06 




42 340,74 

689,00 

1,70 

14 

7Í7lg 

2 pcr u 

-565,52 

978,18 

31,65 




42 342,54 

688,99 

1,70 


A figura 4.11 exibe todas as diferenças dos potenciais das 14 transições que 
contribuem para a formação do perfil de linha de Lyman- 7 , ou seja, AU para cada 
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transiçao. Na figura 4.12 temos todas as transições dos momentos de dipolo elétrico 
para Lyman -7 envolvidas na colisão de prótons com o átomo de H. 



Fignra 4.11: Diferenças dos potenciais das transições em função da distância in- 
teratômica do de Lyman- 7 . 

O perfil de linha calculado para Lyman -7 é mostrado na figura 4.13, para uma 
temperatura de 10 000 K e uma densidade de perturbadores de n p = 10 17 cm -3 . 
Na figura 4.13 pode-se observar muitos satélites no perfil de linha de Lyman 7 , 
entretanto os satélites em À = 998.31 Ã, À = 1 062.64 Ã, A = 688.99 Ã e À = 
733.58 À são os mais proeminentes e, como no caso do satélite em À = 998.31 Ã, 
já observado em espectros de anãs brancas obtidas com o HUT e ORFEUS. 
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Figura 4.12: Transições dos momentos de dipolo elétrico de Lyman -7 do Hij“ em 
fnnção da distância internuclear. O valor de D(R) é dado aqui em unidades atômicas. 
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Figura 4.13: Perfil de linha de Lyman-7 calculado para T=10 000 K e uma densidade 
de prótons de n p — 1, 0 x 10 17 cm -3 . 
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4.3.4 Lyman-(5 

Existem 18 transições moleculares que contribuem para Lyman-á devido a colisões 
do átomo de H com H + . Todas as transições estão listadas na tabela 4.4. 


Tabela 4.4: Transições permitidas e satélites de Lyman-á devido a colisões de H-H + 
para uma distância R entre o par átomo-íon no qual um valor extremo da diferença 
de potencial ocorre. 


N° 

Nível superior 

Nível inferior 

Au h + (cm x ) 

^H+ (Â) 

j?(á) 

1 

5pn u 

lsa g 

- 

- 

- 

2 

6pa u 

ISOg 

- 

- 

- 

3 

QfTTu 

ISOg 

- 

- 

- 

4 

7 / (Tu 

ISOg 

- 

- 

- 

5 

7h7T u 

ISOg 

-679,23 

956,16 

27,56 

6 

8ha u 

lSGg 

- 

- 

- 

7 

00 

e 

lSGg 

-1413,57 

962,93 

36,17 

8 

9 jcr u 

ISOg 

-312,56 

952,82 

44,76 

9 

101 cr u 

ISOg 

-582,87 

955,28 

54,98 

10 

5 SGg 

2 pv u 

- 

- 

- 

11 

6 dog 

2pcr u 

- 

- 

- 

12 

6dn g 

Zpcíu 

- 

- 

- 

13 

7 9 a g 

2 pcr u 

-76,95 

950,69 

26,02 

14 

79*9 

2 pcr u 

- 

- 

- 

15 

8 í(Jg 

2 pcr u 

-1119,44 

960,21 

31,54 

16 

8íl\g 

2 po u 

-99,84 

950,90 

42,38 

17 

9 k(Jg 

2 pc>u 

-1589.40 

964,56 

41,91 

18 

9k7lg 

2 po u 

-475,45 

954,31 

49,00 


A figura 4.14 exibe a diferença da energia potencial de todas as transições que 
contribuem para o perfil de linha de Lyman-á. E na figura 4.15 os momentos de 
dipolo elétrico dessas transições. A figura 4.16 mostra o perfil de Lyman-á calculado 
para a temperatura de T=10 000 K e uma densidade de perturbadores (prótons) 
de riff = 10 17 cm" 3 . O satélite mais proeminente no perfil de linha de Lyman-á 
surge em À = 964.56 Á devido a colisões de átomos de H com prótons formando 
temporariamente a molécula de II.|. 
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Figura 4.14: Diferença da energia potencial das transições permitidas de Lyman-5, 
em função da distância interatômica do H^. 
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Figura 4.15: Momentos de dipolo elétrico das transições de Lyman-á devido ao 
átomo de H perturbado por colisões com prótons. O valor de D(R) está em unidades 
atômicas. 
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Figura 4.16: Perfil de linha de Lyman-5 calculado para T=10 000 K e uma densidade 
de prótons de n p — 1, 0 x 10 17 cm -3 . 
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4.3.5 Série de Lyman 

Como uma síntese dos perfis de Lyman, apresentamos nas figuras 4.17 e 4.18 a soma 
dos perfis de Lya até Lyá. 

Na figura 4.17 são exibidos os perfis que foram calculados, Lya, Ly/3, Lyy e 
Lyô para uma temperatura de T=10 000 K e uma densidade de perturbadores de 
njj = 10 17 cm -3 . Essa figura é similar à figura 3 publicada por Hébrard et al. 
(2003), mostrando que os cálculos e aproximações aqui executadas são consistentes 
com os dados atualmente publicados, mas Hébrard et al. (2003) não havia publicado 
o perfil de linha de Ly<5. A inserção do perfil de linha de Ly ó torna a nossa figura 
mais completa do que a figura 3 de Hébrard et al. (2003). 

A figura 4.18 exibe os perfis teóricos do H^ e os perfis Stark para a temperatura 
T=10 000 K e densidade de perturbadores de n\j = 10 17 cm -3 . Observa-se que o 
perfil gerado pelo alargamento Stark não apresenta os satélites nas asas das linhas, 
que já foram observados. 



Figura 4.17: Perfil total de Lyman-o:, Lyman-/?, Lyman -7 e Lyman-á do H pertur¬ 
bado por prótons ( n^ = 10 17 cm" 3 ) para uma temperatura de T=10 000K. 
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Figura 4.18: Perfil total teórico de Lyman-a, Lyman-/3, Lyman-y e Lyman-5 do H 
perturbado por prótons (n^ = 10 17 cm -3 ) para uma temperatura de T=10 000K do 
IIj juntamente com os respectivos perfis calculados para o alargamento Stark para 
a mesma temperatura e densidade de perturbadores. 
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4.3.6 Balmer-ct 

O perfil total de Balmer-a depende de 32 transições moleculares do íon molecular 
H^, mas nem todas as transições produzem satélites nas asas do perfil de linha de 
Balmer. As transições que geram possíveis satélites nas asas do perfil de linha de Ha 
são listadas na tabela 4.5. Dentre essas transições, as transições 3ga g -4fa u e 4/7t u - 
3dn g são responsáveis pela formação do satélite em À = 8 650 À na asa vermelha do 
perfil de linha de Balmer-a, confirmado via observações em laboratório, ao medir o 
espectro do plasma produzido por laser (Kielkopf et ah, 2002). Na figura 4.19 estão 
plotados os momentos de dipolos das transições moleculares listadas na tabela 4.5 
que contribuem para a formação do perfil de linha de Ha . 


Tabela 4.5: Satélites de Balmer-a devido à colisões de H com H + na distância 
internuclear R para a qual um extremo na diferença dos potenciais da transição 
ocorre,_ 


N° 

Nível superior 

Nível inferior 

Áv H + {cm v 

1 ^ h+ {Ã-) 

R(Ã) 

1 

5 ga-g 

4 fcr u 

-3 844,18 

8 783,59 

13,10 




4956,58 

4954,01 

5,60 

2 

5 gcr g 

2pn u 

-4996,57 

9 772,81 

12,66 

3 

6ha u 

3 d(7g 

40 925,21 

1780,80 

3,17 




-702,54 

6 883,97 

23,03 

4 

4 f-Jiu 

3 d(7g 

-1271,27 

7164,46 

13,13 




2 379,73 

2 562,34 

3,06 

5 

4(1(7 g 

3 pcr u 

-3 214,16 

8 323,01 

6,47 




12145,77 

3 652,99 

1,01 

6 


3d7Tg 

-3 658,85 

8 642,89 

9,33 




7163,74 

4465,72 

1,69 

7 

4/7T„ 

2 SOg 

-5 245,35 

10 016,33 

8,69 




37537,98 

1895,12 

1,00 

8 

5 fcr u 

2 S (7 g 

-424,84 

6 754,83 

15,47 

9 

6 ha u 

2 S (7 g 

-1500,75 

7284,22 

20,33 

10 

4(1(7 g 

2 p-Ku 

-634,82 

6 852,02 

10,20 

11 


2 pn u 

-601.70 

6 836,51 

18,74 

12 

5g°g 

3 po u 

-6 042,60 

10 885,60 

12,08 




23 287,20 

2 596,30 

1,18 

13 

5.9 7T 9 

3pa u 

-1130,91 

7093,14 

16,58 




23 303,15 

2 595,23 

1,18 

14 

6 ha u 

3dra g 

-1117,22 

7086,25 

21,38 




23163,48 

2 604,67 

2,10 
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Figura 4.19: Momentos de dipolo elétrico das transições que contribuem para 
formação dos satélites no perfil de linha de Balmer a, devido ao átomo de Hidrogênio 
perturbado por prótons. As transições plotadas aqui estão listadas na tabela 4.5. O 
valor de D(R) está dado em unidades atômicas. 

A figura 4.20 exibe uma comparação entre o perfil calculado por Koester (2011) 
a partir de uma aproximação de alargamento Stark usando os dados obtidos por 
VCS (Vidal et ah, 1973), para uma temperatura T=12 000 K e densidade de per¬ 
turbadores n^=10 16 cm -3 , e o perfil calculado a partir de uma aproximação quase 
estática, usando potenciais teóricos precisos do íon molecular H^, calculado para 
a mesma temperatura e densidade. Entretanto, há a necessidade de observar que 
a aproximação de alargamento Stark considera a interação envolvendo prótons da 
mesma forma que elétrons e assume um fator multiplicativo de 2 vezes para incluir a 
interação com prótons. Logo, os perfis possuirão uma discrepância da ordem de 0.3 
na escala das ordenadas (logaritmo de 2 na base 10). Observa-se também que em al¬ 
tas energias ambos os perfis coincidem assintoticamente, enquanto, devido a possível 
formação de moléculas, em baixa energia surgem discrepâncias entre os perfis. 

O perfil de linha de Balmer-a foi calculado para uma temperatura de 10 000 K 
e densidade de prótons de 10 1 ' cm -3 e está representado na figura 4.21. 
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MÀ) 


Figura 4.20: Comparação entre o perfil de Balmer-a do 1com o perfil de VCS 
- alargamento eletrônico Stark (Koester (2011)), para T=12 000 K e densidade de 
10 16 cm -3 . 
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Figura 4.21: Perfil de linha de Balmer-a calculado para uma temperatura de 
10 000 K e densidade de prótons de 10 17 cm -3 . Estão indicados na imagem alguns 
dos satélites calculados bem como o satélite A = 8650 Â encontrado em laboratório 
(Kielkopf et ah, 2002) 





Perfis e satélites 


83 


4.3.7 Balmer-/? 

O perfil de linha de Balmer-/3 depende de 46 transições moleculares permitidas do 
H 2 ■ Dentre as 46 transições, 17 transições possuem valores extremos, ou seja, pro¬ 
duzem possíveis satélites no perfil de linha e estão listadas na tabela 4.6, juntamente 
com a distância internuclear onde o satélite é gerado e o respectivo comprimento 
de onda de cada um deles. A figura 4.22 exibe os momentos de clipolo elétrico de 
cada uma das transições listadas na tabela 4.6. É importante salientar que todos os 
valores listados na tabela 4.6 e na figura 4.22 são inéditos, jamais publicados para 
tais transições. Esses valores dos potenciais e momentos de clipolo elétrico do íon 
molecular somente foram publicados em nosso artigo (Santos & Kepler, 2012). 
O perfil de linha calculado de Balmer-/3 é mostrado na figura 4.23 para a temper¬ 
atura T=10 000 K e densidade de perturbadores n p = 1, 0 x 10 17 cm -3 . A figura 4.24 
exibe o perfil de linha de Balmer-/3 para a temperatura T=13 526 K e densidade 
de perturbadores n p = 5, 76 x 10 16 cm -3 de um experimento Z-pinch realizado nos 
laboratórios Sandia (Sandia National Laboratories - USA) como descrito por Falcon 
et al. (2010). A figura 4.25 mostra o perfil de linha de Balmer /3 que calculamos 
para a mesma temperatura e densidade do experimento Z-pinch. Observa-se que o 
perfil teórico possui uma depressão central assimétrica similar ao perfil experimen¬ 
tal mostrado na figura 4.24. De acordo com Wiese et al. (1972), o VCS, ou seja, o 
alargamento Stark não consegue reproduzir a assimetria do pico central, já obser¬ 
vado naquela época, mas somente a depressão central, o que pode ser observado nas 
figuras 4.26 e 4.27. A figura 4.28 exibe o perfil de linha linear de Balmer-/3 teórico 
para a temperatura T=13 526 K e densidade de perturbadores n v = 5, 76 x 10 16 cm -3 
do H 2 somado ao perfil Stark e convoluíclos com o efeito Doppler. A figura 4.29 
mostra a sobreposição das figuras 4.24 e 4.28, onde observa-se que a assimetria no 
experimento Z-pinch é produzida pela formação de quase moléculas 77^ durante a 
colisão entre H-H + . 
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Tabela 4.6: Satélites de Balmer (3 devidos às colisões do átomo de H com H + na 
distância internuclcar R._ 


N° 

Nível superior 

Nível inferior 

Av h + (cm 

^H+ (Â) 

R( A) 

1 

5fn u 

2 S<7g 

-992,00 

5110,59 

13,83 

2 

6hn u 

2 S<7g 

-2 674,11 

5 591,24 

19,90 

3 

5gS g 

2pn u 

-1466,46 

5 237,59 

16,80 

4 


2pvr n 

-1483,16 

5 242,17 

17,74 

5 

6 gcr 9 

3 pcr u 

-2 040,74 

5400,01 

16,50 




23 365,48 

2 276,62 

1,19 

6 

7ia g 

3 pcr u 

-2 977,17 

5 687,62 

24,59 




26 648,74 

2118,28 

1,19 

7 

7m g 

3 por u 

-816,65 

5 065,20 

29,99 




26 649,62 

2118,24 

1,19 

8 

6hn u 

3 d(Jg 

-1811,98 

5 334,12 

21,66 




35 621,24 

1 779,97 

3,17 

9 

8 j(T u 

3da g 

-381,85 

4 956,05 

29,99 




41130,19 

1621,02 

3,19 

10 

5fn u 

3 ditg 

-285,48 

4 932,49 

16,46 




12 088,84 
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Figura 4.22: Momentos de dipolo elétrico das transições responsáveis pela formação 
de satélites no perfil de linha de Balmer-/? devido à perturbação do átomo de H por 
prótons. As transições plotadas aqui estão listadas na tabela 4.6. E os valores de 
D(R) são dados em unidades atômicas. 
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Figura 4.23: Perfil de linha de Balmer-/3 calculado para a temperatura T=10 000 K 
e densidade de perturbadores n p — 1,0 x 10 17 cm -3 . 
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Figura 4.24: Perfil de linha de Balmer-/3 para a temperatura T=13 526 K e densidade 
de perturbadores n p = 5, 76 x 10 16 cm” 3 de um experimento Z-pinch, como observado 
por Falcon et al. (2010). 
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Figura 4.25: Perfil de linha de Balmer-/3 teórico somente do efeito de H — H + para 
a temperatura T=13 526 K e densidade de perturbadores n p = 5, 76 x 10 16 cm -3 . A 
janela no topo direito mostra uma expansão do pico central, mostrando a assimetria. 
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Figura 4.26: Perfil de linha linear de Balmer-/3 teórico para a temperatura 
T=13 526 K e densidade de perturbadores n p = 5, 76 x 10 16 cm -3 do H^ juntamente 
com o perfil Stark. 
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Figura 4.27: Perfil de linha de Balmer-/3 teórico para a temperatura T=13 526 K e 
densidade de perturbadores n p = 5, 76 x 10 16 cm -3 do juntamente com o perfil 
Stark. 


















F(cm 


Perfis e satélites 


91 


0,02 


0,016 


0,008 


0,004 


Q I _I_I_I_I_I_I_ 

4750 4800 4850 4900 4950 

X(À) 

Figura 4.28: Perfil de linha linear de Balmer-/3 teórico para a temperatura 
T=13 526 K e densidade de perturbadores n v = 5, 76 x 10 16 cm -3 do somado 
com o perfil Stark convoluídos com o efeito Doppler. 
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Figura 4.29: Perfil de linha de Balmer-/3 para a temperatura T=13 526 K e densidade 
de perturbadores n p = 5, 76 x 10 16 cm -3 de um experimento Z-pinch, como observado 
por Falcon et al. (2010) e sobreposto ao mesmo o perfil de linha linear de Balmer-/3 
teórico para a temperatura T=13 526 K e densidade de perturbadores n p = 5, 76 x 
10 16 cm -3 do H 2 juntamente com o perfil Stark convoluídos com o efeito Doppler. 
As linhas Stehlc e Finer Grid foram calculadas por Michael ff. Montgomery com o 
programa hlinfitJemke.f90 . 
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4.3.8 Balmer -7 

O perfil total da linha de Balmer -7 depende de 60 transições moleculares do ■ 
O perfil de linha teórico mostrado na figura 4.30 foi calculado para a temperatura 
T=10 000 K e para a densidade de perturbadores (prótons) n p = 1,0 x 10 17 cm -3 . 





Figura 4.30: Perfil de linha de Balmer -7 calculado para a temperatura T=10 000 K 
e para uma densidade de perturbadores n p — 1,0 x 10 17 cm -3 . 

Observando a figura 4.31, a assimetria é clara e vai alterar os cálculos de opaci¬ 
dades. 
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Figura 4.31: Perfil de linha linear de Balmer -7 calculado para a temperatura 
T=10 000 K e para uma densidade de perturbadores n p = 1,0 x 10 17 cm -3 do H^ 
em comparação com o perfil Stark para mesma temperatura e densidade de pertur¬ 
badores. 
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4.3.9 Balmer-d 

O perfil de Balmer-á é obtido a partir de 74 transições moleculares do íon H^- O 
perfil de linha teórico foi calculado para a temperatura de 10 000 K e uma densidade 
de rip = 1,0 x 10 17 cm -3 perturbadores (prótons) e está exibido na figura 4.32. 



Figura 4.32: Perfil de linha de Balmer-á calculado para a temperatura T=10 000 K 
e para uma densidade de perturbadores n p — 1,0 x 10 17 cm -3 . 
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4.3.10 Balmer-e 

O perfil de linha total de Balmer-e depende de 88 transições moleculares do íon 
■ A figura 4.33 exibe o perfil teórico calculado para a temperatura T=10 000 K 
e densidade de perturbadores n p = 1.0 x 10 17 cm -3 . 



X(À) 


Figura 4.33: Perfil de linha de Balmer-e calculado para a temperatura T=10 000 K 
e para uma densidade de perturbadores n p = 1.0 x 10 17 cm -3 . 
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4.3.11 Balmer-8 

O perfil de linha total de Balmer-8 é gerado a partir de 102 transições moleculares 
do H 2 ■ Seu perfil teórico é mostrado na figura 4.34 e o mesmo foi calculado para 
uma temperatura T=10 000 K e densidade de perturbadores n p = 1,0 x 10 17 cm -3 . 



X(À) 


Figura 4.34: Perfil de linha de Balmer-8 calculado para a temperatura T=10 000 K 
e para uma densidade de perturbadores n p — 1,0 x 10 17 cm -3 . 
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4.3.12 Balmer-9 

O perfil de linha total de Balmer-9 depende de 116 transições moleculares do íon . 
Seu perfil teórico foi calculado para a temperatura de T=10 000 K e densidade de 
perturbadores n p = 1.0 x 10 17 cm -3 . A figura 4.35 exibe o perfil teórico de Balmer-9 
obtido. 





Figura 4.35: Perfil de linha de Balmer-9 calculado para a temperatura T=10 000 K 
e para uma densidade de perturbadores n p = 1.0 x 10 17 cm -3 . 
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4.3.13 Balmer-10 

O perfil de linha total de Balmer-10 é obtido a partir de 130 transições moleculares do 
íon ■ O perfil teórico é plotado na figura 4.36 e foi calculado para a temperatura 
de T=10 000 K e densidade de perturbadores n p = 1,0 x 10 17 cm -3 . 



X(À) 


Figura 4.36: Perfil de linha de Balmer-10 calculado para a temperatura T=10 000 K 
e para uma densidade de perturbadores n p = 1.0 x 10 17 cm -3 . 
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4.3.14 Série de Balmer 

Como uma síntese dos perfis de Balmer, apresentamos nas figuras seguintes a soma 
dos perfis de Ha até H10, perturbados por colisões de hidrogênio com prótons para 
a temperatura de T=10 000 K e número de perturbadores n p = 1,0 x 10 17 cm -3 . 

A figura 4.37 exibe os perfis lineares dessa série. 



Figura 4.37: Contribuição do H — H + para o perfil de linha linear total de Balmer 
calculado para a temperatura T=10 000 K e para uma densidade de perturbadores 
n p = 1 , 0 x 10 17 cm -3 . 

Na figura 4.38 colocamos os perfis Stark juntamente com os perfis obtidos para 
oH+, para comparação. 

A figura 4.39 exibe os perfil de Balmer-a, Balmer-/?, Balmer- 7 , Balmer-Ú, Balmer- 
e, Balmer- 8 , Balmer-9 e Balmer-10 juntos para a mesma temperatura, ou seja, 
T=10 000 K e a mesma densidade de perturbadores np = 10 17 cm -3 . A normal¬ 
ização do perfil da figura 4.39 foi calculado para o perfil somado e não individual 
como na figura 4.38. 
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Figura 4.38: Contribuição do H — H + e do elétron (Stark) para o perfil de linha 
total de Balmer calculado para a temperatura T=10 000 K e para uma densidade 
de perturbadores n p — 1, 0 x 10 17 cm -3 juntamente com o perfil de linha Stark, para 
a mesma temperatura e densidade de perturbadores (elétrons para o Stark). A 
normalização foi calculada para cada transição individual da série de Balmer. 
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Figura 4.39: Contribuição do H—H + para o perfil de linha total de Balmer calculado 
para a temperatura T=10 000 K e para uma densidade de perturbadores n p — 1, 0 x 
10 17 cm~ 3 . A normalização foi calculada para a soma das transições da série de 
Balmer. 





Capítulo 5 

Discussão e Conclusão 


Teorias satisfatórias e dados para o perfil de linha existiam somente para o alarga¬ 
mento Stark do hidrogênio neutro (Lemke, 1997, Vidal et ah, 1973) e para 21 linhas 
óticas do Hélio neutro (Barnard et ah, 1969, Beauchamp et ah, 1997). Somente 
as três primeiras linhas de Lyman do átomo de hidrogênio alargadas por pertur¬ 
badores neutros e ionizados, incluindo alguns satélites característicos, eram bem 
descritas pelo trabalho de Allard et al. (2009) e incluidos nos cálculos de modelos 
atmosféricos. Os dados atômicos podiam ser obtidos somente a partir de bancos de 
dados, principalmente das listas de linhas de (Kurucz & Bell, 1995) e do banco de da¬ 
dos VALD ( Vienna Atomic Line Database ) Kupka et al. (1999, 2000), Ryabchikova 
et al. (1997), Piskunov et al. (1995). Entretanto, para os outros casos a situação 
não era satisfatória, de modo que havia uma grande necessidade de cálculos e dados 
que possam explicar de maneira adequada os perfis de linha e as linhas satélites 
observáveis. 

Devido a essa inexistência, calculei os potenciais e momentos de dipolo elétrico 
do íon molecular , gerado por colisões entre átomos de hidrogênio e prótons, até 
o nível n=10 (número quântico principal) do mesmo. Calculei também os perfis 
de linha das séries de Lyman e Balmer com densidades apropriadas para aplicação 
em atmosferas de estrelas anãs brancas [tipicamente n = 10 17 cm -3 , Allard et al. 
(1994)], dados e cálculos nunca publicados até nosso artigo (Santos & Keplcr, 2012), 
e disponibilizados em: 

http: / / astro.if.ufrgs.br / márcios. 

As observações em laboratório de plasma produzido por laser confirmaram que a 
maioria dos satélites que aparecem em Lyman-a são devidos às colisões entre átomos 
de hidrogênio neutros e prótons. Essa confirmação experimental da teoria reforça a 
identificação desses satélites em espectros de Lyman-o: de estrelas anãs brancas e de 
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estrelas A Bootis. Nos plasmas em altas temperaturas, tais como plasmas produzidos 
por laser, processos de recombinação dieletrônicos são muito importantes. Esse 
processo afeta fortemente a abundância de íons e produz linhas satélites dieletrônicas 
mas as variações temporais são de difícil cálculo. As linhas satélites de emissão a 
partir de íons altamente ionizados são frequentemente utilizadas para diagnósticos 
de diversos tipos de plasmas, por exemplo, TOKAMAK , fusão inercial, astrofísico 
etc. 

A comparação das asas de perfil de Lyman-a com perüs de linha gerados por 
modelos é uma ferramenta para determinar as densidades de átomos neutros e ioniza¬ 
dos em plasmas hidrogênicos. Esses experimentos também confirmam que a variação 
do momento de dipolo radioativo é um fator importante na determinação das asas 
distantes do perfil de Lyman-a; quando o momento de dipolo D(R) difere signi¬ 
ficativamente do seu valor assintótico para uma distância R próxima à região de 
formação de um satélite, a intensidade da asa pode ser aumentada ou diminuída 
consideravelmente. 

Então, num espectro de absorção induzida por colisão (do inglês Collision In- 
duced Absorption - CIA), tal como o mencionado, o momento de dipolo pode ser 
mais importante do que o potencial na determinação da forma de um satélite. Tais 
perturbações de muitos corpos alteram a linha da asa distante do perfil, adicio¬ 
nando vários satélites e produzindo um contínuo forte a partir do ultravioleta para 
o visível. Consequentemente, o estudo de atmosferas de estrelas anãs brancas, plas¬ 
mas e experimentos associados (Z-pinch) é fundamental para a interpretação dos 
dados observacionais e representa a ligação dos modelos teóricos com a análise de 
seus resultados, ou seja, a criação de modelos adequados para explicar as observações 
de tais processos físicos e astrofísicos. 

Os espectros observados de estrelas anãs brancas exibem uma forma carac¬ 
terística complexa na qual intervem processos de difusão, acreção e convecção. Até 
o momento, essa forma característica não é compreendida em todos os seus detalhes, 
devido a inexistência de um modelo atmosférico que inclua todos os processos físicos 
relevantes que ocorrem na atmosfera. Seu estudo é a chave na interpretação de es¬ 
pectros e na análise das predições da teoria de evolução. Portanto, o assunto desse 
trabalho também é de real interesse para modelamento dos espectros observados. 

A aproximação teórica utilizada foi baseada na absorção quase molecular do 
devido a interação H — II ' , calculada na aproximação quase estática. Nessa aprox¬ 
imação, além de fatores estatísticos, os perüs de linha dependem principalmente 
da derivada dVij/dR da diferença de energia entre os estados moleculares i e j de 
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diferente paridade e dos momentos de dipolo moleculares correspondentes. 

O presente trabalho foi divido essencialmente em duas grandes partes: ( 1 ) os 
cálculos dos dados moleculares e ( 2 ) o cálculo dos perfis correspondentes. 

Na primeira parte, determinamos os potenciais e clipolos da quase molécula 
pois, apesar de já ter sido estudada por diversos autores, não haviam dados molecu¬ 
lares para estados eletrônicos superiores a n = 3 disponíveis na literatura. Apresen¬ 
tamos o método usado para obter potenciais e momentos de dipolo elétrico precisos, 
bem como os dados moleculares para estados até n=10. Os potenciais e momentos 
de dipolo elétrico obtidos foram comparados com os existentes publicados até o pre¬ 
sente momento, por Madsen & Peek (1971), Ramaker & Peek (1972) e Ramaker & 
Peek (1973), coincidindo com excelente precisão. 

A segunda parte foi dedicada aos cálculos dos perfis de linhas das séries de 
Lyman e Balmer, onde usei a aproximação quase estática (Bates, 1952, Rohrmann 
et ah, 2011). Nessa aproximação, os valores extremos da derivada dVij/dR produzem 
divergências no perfil de linha, e é isso que leva ao surgimento de satélites nas asas 
dos mesmos. 

Ao longo dos cálculos, foi mantida a máxima resolução possível para não perder 
a forma global correta das asas das linhas e nem as posições exatas teóricas dos 
satélites ao longo do perfil. Calculei a posição de todos os satélites das principais 
transições das séries de Lyman e Balmer, o que pode ser verificado tanto nos perfis 
de linha obtidos quanto especificamente nas tabelas 4.1, 4.2, 4.3, 4.4, 4.5 e 4.6. 

Na figura 4.5 apresentamos individualmente cada uma das 6 transições molec¬ 
ulares do H 2 que contribuem para a formação do perfil de linha de Lyman-o:, e 
pode-se perceber que a transição entre os níveis 2 sa g e 2 pa u gera uma asa vermelha 
distante que se estende até a parte azul do espectro ótico, que pode ser a fonte de 
opacidade que faltava nos modelos de atmosferas de estrelas anãs brancas, corrobo¬ 
rando com a ideia de Kowalski & Saumon (2006). Os perfis de Lyman-/? e Lyman -7 
estão sendo estudados desde Allard & Kielkopf (1982), entretanto somente em Al- 
lard et al. (2009) eles incluiram os cálculos do momento de dipolo em Lyman- 7 , 
mas que não foram publicados. O perfil de Lyá teórico que obtive nunca havia sido 
calculado até então, e apresenta um satélite proeminente em À = 964.56 Ã. 

O perfil de linha de Balmer-o: é mostrado na figura 4.21, com melhor resolução 
que os artigos previamente publicados. Para o perfil de linha de Balmer-/?, encontrei 
os possíveis satélites listados na 4.6, sendo que alguns podem ser verificados na 
figura 4.23 para a temperatura de 10 000 K e uma densidade de prótons de n v = 
1.0 x 10 17 cm -3 . Observa-se também que esse perfil possui uma depressão central 
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assimétrica similar ao perfil experimental da figura 4.24 obtido por Falcon et al. 
(2010) junto aos Laboratórios Sandia. Sabe-se que o alargamento Stark (Wiese 
et al., 1972) não reproduz a assimetria na altura dos picos, mas é causado pelos 
prótons em colisão, formando a quase molécula de ■ Os demais perfis de linha 
da série de Balmer foram apresentados nas figuras de 4.30 até 4.36, para a mesma 
temperatura e densidade de perturbadores. 

Esse trabalho mostrou que os perfis obtidos por colisões com íons diferem sig¬ 
nificativamente do alargamento Stark por elétrons, mostrando não apenas linhas 
satélites, mas também valores significativos das asas dos perfis que afetam a opaci¬ 
dade. Está centrado no estudo teórico dos perfis de linha do hidrogênio perturbado 
por colisões com prótons, representando um avanço significativo no campo da mod¬ 
elagem de perfis e de atmosferas estelares de anãs brancas DA, gerando perspectivas 
para novos estudos e pesquisas da etapa final de evolução da maioria das estrelas, 
que pode revelar maiores informações sobre a história da nossa galáxia. 

A continuação desta pesquisa envolverá a comparação com perfis de densidade 
única medidos em Sandia, inclusão destes perfis no cálculo de modelos de atmosfera, 
bem como o cálculo similar para o H 2 . 
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6.1 Unidades atômicas 

Para os cálculos deste trabalho, foi conveniente fazer uma troca de unidades para 
unidades atômicas (u.a.), também chamadas de unidades de Hartree: a unidade de 
massa é a massa do elétron e a unidade de carga é a carga do elétron, de momento 
angular é h. 

A unidade atômica de comprimento é 

bohr = a 0 = 0, 52918Â. 

Quando as unidades foram mudadas para unidades atômicas, h, m e e, trocamos o 
valor dessas constantes para o valor numérico 1. Assim, a utilização destas unidades 
simplifica também o hamiltoniano: 

H 2 2 1 e 2 1 2 1 

2 m 47r6 o r 2 r 

e a unidade atômica de energia recebe o nome de hartree: 

e 2 

— = 27, 212eV = l hartree 
ao 
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Símbolo 

Quantidade 

S.I. 

u.a. 

m 

massa (do elétron) 

9,1094xl0~ 31 kg 

1 

e 

carga (do elétron) 

1,6022x10” 19 C 

1 

h 

Constante de Plank (h/2n) 

6,626 x10” 34 Js/2tt 

1 

o 0 

raio de Bohr 

0,52918x 10 _1 °m 

1 

4vre 0 

Permissividade do vácno 

l,113x 10” 10 C 2 /Jm 

1 

c 

velocidade da luz 

2,998 x 10 8 m/s 

- 

k 

Constante de Boltzmann 

l,38066x 10 _23 J/K 

- 

N a 

número de Avogadro 

6,02205 x 10 23 mol _1 

- 
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6.2 Artigo publicado 

Theoretical study of the line profiles of the hydrogen perturbed by colli- 
sions with protons 

Monthly Notices of the Royal Astronomical Society, Volume 423, Issue 1, 
pp. 68-79. Publicado em Jun 2012 DOI: 10.1111/j.1365-2966.2012.20631.x 

We present theoretical calculations of the quasi-molecular line profiles for the Ly- 
man ( Lya, Lyf3, Ly 7 and Lyô ) and Balmer ( Ha , Hf3, H 7 , H5, He, H 8 , H 9 and 
H 10) series pertnrbed by collisions with protons. In all calculations, we include the 
dependence of the dipolc moments as a function of the internuclear distance clur- 
ing the collision. The broadening from ion collisions must be added to the normal 
clectron Stark broadening. 


Keywords: Astronomy: line: profiles, white dwarfs 
arXiv: arXiv:1204.4769 

http://onlinelibrary.wiley. com/doi/10.1111/j. 1365-2966.2012.20631.x/full 
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